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ABSTRACT 


Interoperability of commercial Land Mobile Radios (LMR) and the militarw’s tac- 
tical LMR is highly desirable if the U.S. government is to respond effectively in a na- 
tional emergency or in a joint military operation. This ability to talk securely and 
immediately across agency and military service boundaries is often overlooked. One way 
to ensure interoperability is to develop and promote federal communications standards 
(FS). 

This thesis surveys one area of the proposed FS 1024 for LMRs; namely, the error 
detection and correction (EDAC) of the message indicator (MI) bits used for 
cryptographic synchronization. Several EDAC codes are examined (Hamming, Quad- 
ratic Residue, hard decision Golay and soft decision Golav), tested on three FORTRAN 
programmed channel simulations (INMARSAT, Gaussian and constant burst width), 
compared and analvzed (based on bit error rates and percent of error-free superframe 
runs) so that a “best” code can be recommended. Out of the four codes under study, the 
soft decision Golay code (24.12) is evaluated to be the best. This finding is based on the 
code’s ability to detect and correct errors as Well as the relative ease of implementation 
of the algorithm. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. BACKGROUND 
1. Why Land Mobile Radio Standards are Needed 


Interoperaivility, connectivity and security of the U.S. land mobile radio 





(LMR) communications systems are important objectives at all levels - federal, 
State, and local. In many ways the U‘S. falls short in meeting these communi- 
cations objectives. “This situation is even more serious when viewed in the Na- 
tional Security Emergency Preparedness (NSEP) context’[Ref. 1]. During such 
an emergency, private users of a land mobile public safety svstem may need to 
communicate with federal! defense and law enforcement agencies to effectively 
execute emergency operations. This is not currently possible. 

Within the defense community, interoperability between various services 
is a longstanding problem. An example that illustrates the extent of this problem 
is provided by the communication problems encountered during the Grenada 
conflict. Many rumors have circulated about a soldier who used his AT&T card 
to call back to the continental United States by phone, since he could not com- 
. municate with soldiers on the other side of the island due to lack of radio inter- 
operability.! For all future joint military conflicts, interoperability between the 
services IS a necessity. 

The absence of a federal communications standard for the LMR, until 
most recently, has encouraged the proliferation of a diverse group of technically 
incompatible LMRs. The fact that the commercial market h’s produced a wide 
varicty of state-of-the-art equipment is good in one sense -- there is much to 
choose from and th. competition in the commercial market place keeps the cost 
low. However, each manufacturer produces his own “best” product with no 
standards or guidelines to follow. Modulation schemes, for example, may be 


different. This makes communication possible only with someone who has pur- 


1 This story has never been substantiated. 











chased a radio from the same vendor. But if the overall objective is to promote 
LMR interoperability across agency and service boundaries, then more thorough 
Standardization is necessary. 

2. Current and Proposed LMR Standards. 

There is presently a LMR standard in place -- Federal Standard (FS) 
1023. Because of it’s late adoption (September 1989), it did not have an impact 
on existing defense or civilian LMR systems; many incompatible LMRs were al- 
readv fielded. FS 1023’s near-term objective was to prevent further proliferation 
of incompatible systems until a new standard for future upgraded systems could 
be agreed upen. 

FS 1023 will be obsolete in the near future, in part because ”...increasing 
demands on the fixed spectrum availability along with overseas evolution towards 
nairower channels have caused the need for efficient spectrum use.” [Ref. 2] 
Today. new technological advancements enable more efficient use of the fre- 
quency spectrum. FS 1023 is based on 25 kHz allocated channel spacing: the 
proposed standard, FS 1024, will require a narrower allocation of frequency 
channel spacing -- 12.5 or 6.25 kHz. 

LMR interoperability, through the adoption of FS 1024. is an achievable 
mid-term goal. Such an LMR standard would help to alleviate both interoper- 
ability and spectral congestion concerns. 

Because of the interest and justifed need for a future federal standard, the 
National Communications System (NCS) office of Technology and Standards, 
the National Security Agency (NSA) and the Land Mobile Radio (LMR) sub- 
committee of the Federal Telecommunication Standards Committee (F TSC) are 
working together to develop Federal Standard 1024 to meet the future needs of 
all LMR users -- federal (civil and defense), state, and local. The official draft 
FS !024 is expected to be released in Spring 1990 for 90 day industry coordi- 
nation and comment [Ref. 3]. Other correspondents involved include the US. 
Department of Justice (DOJ), the Assistant Secretary of Defense for C31, and the 
Gencral Electric Company [Refs. 1.4, 5. 6. 7, 8. and 9]. 





3. Advantages and Disadvantages of FS 1024 

There are three clear advantages with the adoption of FS 1024. First, full 
and open competition between commercial LMR vendors will become healthier. 
Clearly the company with the best product for the money will win a competitive 
procurement. Second, a healthy market will produce a greater availability of 
LMR products. Third, it is likely that the standard will extend beyond the 
bounds of the intended commercial LMR market to the tactical defense radio 
market. This “status quo” effect, driven by the adoption of FS 1024, would have 
a positive influence on intra-service communications compatibility and interop- 
erability. 

Most often, each service provides the communications for its own units. 
There is a wide variety of tactical radios fielded for dozens of unique applications. 
It is more of an exception when these radios are able to interoperate. FS 1024 
may be seen, then, as a first step in enhancing cross-service interoperability re- 
quirements for tactical (land mobile) units. 

The main disadvantage of FS 1024 is that standardized radios would not 
be backwards compatible with existing radios that comply with FS 1023. An 
upgrading of all existing systems, a high-cost and highly unlikely proposition. 
would be necessary before full interoperabilitv could be obtained. Over time this 
disadvantage will be lessened since an in-place Federal Standard for new LMR 


systems Will ultimately drive the replacement of existing systems. 


B. DESCRIPTION OF LAND MOBILE RADIO FS 1024 
The narrowband digital LMR standard includes three criteria to ensure 
interoperability. These three criteria are: 
® voice coding (digitization) technique. 
e the radio freyuency (RF) modulation technique, 


e the cryptographic algorithm and cryptographic synchronization. 














The voice coding technige refers to another federal standard - FS 1016.2 This 
standard specifies CELP (codebook excited linear predictive) coding as the re- 
quirement. CELP coding will be implemented with a 8 KHz sampling rate. The 
algorithm itself contains three functions: short term spectral prediction, long de- 
lay adaptive codebook “pitch” searches, and innovative stochastic codebook 
search. 

The second criterion, RF modulation, is still under consideration. Possibil- 
ities include 4-ary frequency shift keying (FSK), tamed frequency modulation 
(FM), quadrature differential phase shift keying (QDPSK) and 2-4 shift 
QDPSk.3 Final selection will be based on spectral efficiency (minimizing adjacent 


channel interference) and power efficiency. 


C. SCOPE AND GOAL OF THESIS STUDY 

The third criterion, cryptographic algorithm and sychronization is the portion 
of FS 1024 on which this thesis will concentrate. The focus is on the error de- 
tection and correction (EDAC) scheme for the message indicator (M]) portion of 
the synchronization bits. within the allocated transmission format. 

The goal of this study is to determine the most suitable EDAC code, among 
four -- Hamming. Quadratic Residue Code (QRC), Golay hard decision and 
Golay soft decision -- for cryptographically securable LMRs. Factors affecting 
the findings are also discussed. Analyses and findings are based on a computer 
simulation. The thesis only addresses block codes. Convolutional codes Were 
eliminated from consideration by DIRNSA based on a study by NASA [Ref. 
11]. Longer processing delays -- characteristic of convolutional codes -- are con- 
sidered unacceptable for use in LMIR, where a more immediate response is re- 
quired. 


2 FS 1016 is proceeding through the FISC approval process [Ref. 10: p. 3] and is expected to 
be approved by the Summer of 1990. 


3 In an analysis performed by NSA. coherent and mainly noncoherent x.4 shift QDPSK were 
used. [Ref. 2] 





D. ORGANIZATION OF THESIS 

Chapter II discusses the transmission frame (bit) structure and proposed 
interleaving scheme. Interleaving and EDAC serve to mitigate noise induced er- 
rors. 

Chapter II] discusses and describes the EDAC codes considered. 

Chapter IV describes the simulation and tests performed. 

Chapter V presents the test results. 

Chapter VI provides conclusions and recommendations for further study. 

The appendix contains the FORTRAN code for the computer simulations 
used (Appendix A), the selected output data of the simulation runs (Appendix 
B). and the codeword interleave tables for each code type (Appendix C). Golav 
hard and soft decision use the same interleave table, therefore, there are only 


three algorithm ‘types’ out of the four codes tested. 


E. THESIS COLLABORATION AND SUPPORT 

This thesis was performed in consonance with DIRNSA R556 as an aid in 
their study and analvsis of the FS 1024 issues for DIRNSA, V2 (secure voice ac- 
quisition) and the National Communications System (NCS). Also, GTE Gov- 
ernment Systems Corporation - Electronic Defense Communications Division in 
Waltham. Massachusetts is currently evaluating portions of the new FS 1024 for 
DIRNSA under contract number MSA904-90-C-6014. 


a 


II. TRANSMISSION FRAME STRUCTURE AND INTERLEAVING 


A. FS 1024 SUPERFRAME AND FRAME 
Figure | on page 7 shows the bit allocation for the FS 1024 data rate of 8000 

bits:s. Table 1 on page 8 summarizes bit allocation and the respective bit rate 
for each superframe. For each second, 2400 are Error Detection and Correction 
(EDAC, or “parity”) for the information bits, 4800 of the bits are information bits 
(CELP processed) and the remaining 800 bits are used for overhead processing. 
Bits allocated for overhead processing include bits for: 

® framing 

® mode control 

@ system control 


® cryptographic synchronization (i.e., message indicator (MI) bits) 


OM (800 BITS/S) INCLUDES: FRAMING, CRYPTO SYNC, MODE CONTROL, SYSTEM CONTROL 





Figure 1. Data Rate Allocation (Source: GTE Government Systems Corp., 
LMR program review, Waltham, MA) 





Table 1. SUPERFRAME BIT COUNT (420 MS) 


Voice EDAC (parity) 2400.00 
Voice (CELP bits) 4800.00 
Framing 114.28 

Mode Control (MC) 9.52 
MC EDAC 23.80 

Msg Indicator (MI) 171.42 
MI EDAC 171.42 


Reserved bits 14.29 


System Control 295.24 





Adapted from GTE - Government Systems Corporation 


For the recommended superframe transmission and interleaving format, there 
are 8 modes of operation defined by 4 bits in the 800-bit overhead contro} field: 
four are for the encrypted mode and four for the plain text mode. Only the en- 
crypted voice mode, designated as 0100 in the mode control block, will be ad- 
dressed in this thesis since the message indicator (MI) bits are not used in the 
unencrypted or plain text frame structure. 

Figure 2 on page 9 shows the proposed frame structure for both the en- 
crypted and the plain text mode. The shaded portion, the information super- 
frames, is the only section of the message considered. The plain text mode does 
not incorporate EDAC or MI bits; thus, plain text superframes are half the du- 
ration. 

Figure 3 on page 10 provides a one second snapshot of 8000 bits. Each 
superframe contains fourteen 30 msec frames. A total of 72 MI bits plus 72 MI 


parity bits will be interleaved throughout these 9 of these 14 frames, in a process 


T = 420 ms (ENCRYPTED TEXT MODE) 
T = 240 ms (PLAIN TEXT MODE) 


_TEADER INFORMATION 


Ene AND SUPER - SUPER 


eon ms “ome ms 
O14(5a)-6/21/89-C8D 





Figure 2. Transmission Format (Source: GTE Government Systems Corp., 
LMR program review, Waltham, MA) 


explained below. Three levels are depicted in Figure 3 on page 10. The top level 
depicts 8000 bits or 8000/3360 = 2.38 superframes. The next level shows the 
superframe broken down into the fourteen 30 msec frames. Finally in the lowest 
level, one of these 240 bit frames shows the arrangement of the 16 interleaved bits 
within one frame. Of the 14 frames, frames 4 through 12 (only) contain the MI 
bits and MI EDAC (or “parity”) bits. 


B. PROPOSED INTERLEAVING 

The process of interleaving involves dispersing critical bits over the super- 
frame, in regular fashion, so that there is a smaller chance of these bits being 
corrupted by noise bursts. Interleaving ensures, for example, that the 72 MI bits 
and the 72 MI parity bits are spread in a pre-determined arrangement through- 
out the superframe. In this way a long noise burst will not corrupt these critical 
portions of the message. If the bits were grouped together, even a 10 ms burst (80 
bits) could destroy the usefulness of the entire MI, and prevent successful 


cryptographic synchronization. Instead, by using an interleaving technique, only 
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a few of the critical MI bits may be corrupted; these should be detected and cor- 
rected by the interleaved MI parity bits. 4 

An evaluation was performed by GTE on two proposed interleaving schemes 
(I and II) for the FS 1024 superframe format [Ref. 12]. The testing was per- 
formed using mobile vehicles communicating in an urban environment. Scheme 
Il outperformed Scheme I for vehicle speeds approximately three times slower in 
the five LMR frequency bands (132-172MHz, 406-420MHz, 450-512MHz, 
806-512MHz and 851-866MHz).5 This is because the EDAC codewords (MI bits 
+ parity bits) are more evenly interspersed over each frame in the Scheme Il 
interleaver. The subsequent description and analysis address interleaver scheme 
II only, since FS 1024 will likely incorporated that scheme. 

Figure 3 on page 10 also shows the interleaving scheme used by the EDAC 
computer simulation described in Chapter 1V. Frames 1 through 3 are reserved 
for framing and mode control bits; Frames 13 and 14 contain system control bits. 
Every frame contains the voice - CELP processed information bits and its parity 
bits. The interleaving of MI bits and MI parity bits begins at the beginning of 
frame #4 or bit 721 and alternates skipping first 10 bits then 20 bits then 10 bits. 
etc. all the way through the frame #12 ending at bit 2680. Appendix C includes 
three codeword interleaving tables, one each for the Golay (24.12) code. the QRC 
(48,24) code and for the Hamming (8.4) code. The MI bits plus their parity bits 
(or “codewords”) are broken up differently for the Golay, QRC, and Hamming 
codes to achieve maximum spread over frames 9 thru 12, the interleaving frames. 
For QRC the codewords are broken up into 4 bit chunks. for Golay the 
codewords are broken up into 2 bit chunks and for Hamming | bit chunks. 

Figure 4 on page 13 illustrates frame composition for the QRC (48,24). For 
the QRC there are 3 codewords required (3 x 48 bits = 144 bits). Each codeword 


4 More critical. however, is the interleaving of the MI parity bits. For if non-interleaved panty 
bits were corrupted by a noise burst. then MI bits also corrupted by another random noise burst 
would have no chance of being detected and corrected by the chosen EDAC scheme. At the verv 
least. then. MI parity must be interleaved. Because of the random nature of noise however, both 
MI bits and their parity bits are both. optimally. interleaved. 


3 For both schemes the lowest frequency band had the worst performance. 











block pictured is a piece of one of the 3 codewords taken 4 bits at a time. A!l bits 
other than these 144 out of 3360 bits’superframe are treated as “don’t care” 
Within the analvsis program. The EDAC protection of the information bits 
(CELP processed voice) was analyzed by NSA, R556 using Golay (hard and soft 
decision) code [Ref. 2]. 

In Chapter III, block codes and the codewords used for each code will be 
discussed in more detail. Most importantly, each code will be examined in terms 


of its ability to detect and correct errors. 
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Figure 4. Superframe Composition (Source: GTE Government Systems Corp., 
LMR program review, Waltham, MA) 
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Il. EDAC CODES EXAMINED 


A. INTRODUCTION 

All of the codes examined in this thesis -- Golay(24,12), QRC(48,24), and 
Hamming(8&.4) are either in the class of (n,k) block repetition codes or (n,k) cyclic 
block codes. For n total bits:block and k information bits/block, the number of 
parity bits equals (n-k). and code rate R = kn. All of the codes have a code rate 
R = 1 2, which lends itself to easier bit manipulation. 

For block codes. a generating polynomial, g/x), is used to determine the 
parity checks performed. The number of parity bits determine the the number 
of checks performed. This parity check group is formed into a check matrix or 
H matrix for ease of manipulation. 

The n-bit blocks (or “codewords”), formed by the product of the encoder (the 
parity bits) and information bits, becomes the transmit vector, }’/x). The re- 
ceived vector, before decoding. is designated R/x). 

Prior to decoding. a syrdrome is used as a binary number block by the code 
algorithm to mark errors and the position of the errors for correction. If the 
syndrome equals zero, then there are no errors. The decoding and correction of 
the bits are algorithm dependent. Severa! different decoding schemes may be 
possible for each code.6 

The number of possible errors corrected for a code is determined not only by 
the codeword size (n) and the information stream length (k). but by the minimun 
Hamming distance. dp. This distance is the minimum number of bits that differ 
among all pairs of codewords in the code set [Ref. 16: p. 24]. The code can cor- 


rect | 2(d,,,) number of errors if d.,, is even and 1,2(d,,, - 1) errors if d,,, is odd. 


B. GOLAY CODES 
The perfect cyclic Golay (23,12) code is a triple-error correcting code, and can 


be implemented in software. Typically, one extra bit is added to parity yielding 


6 For derivations and more study on linear block codes refer to [Refs. 13. 14, and 15]. 








the Golay (24,12) code. This is done to maintain code rate = 1.2 and ease of 


algorithmic manipulation. The generating polynomial used is: 
: gxpaxtn xr gx 4 x84 xP txt 


The Golay (24.12) code can be improved by utilizing analog informat on 
asssociated with the channel bits. In this wav the bit error rate in white Gaussian 
noise can be improved about 2 dB. {[Ref. 17] This “soft” decision-inaking uses the 
Value of the previous baud to calculate confidence intervals. These confidence 
values are rea/ values rather than actual baud values. This concept is illustrated 
by the constellation. Figure 5 on page 16 [Ref. 2). 

The X and Y axes of Figure 5 on page 16 represent the least significant bit 
(Isb) and the most sign: icant bit (msb) decision boundaries of the two baud being 
compared. The dot represents the unquantized phase difference between present 
baud and previous baud and the circle represents the magnitude at the present 
baud relative to the low-pass filter receive baud magnitude. The confidence for 
each sb is then the phase difference from the dot to its corresponding nearest 
boundary. Fades are flagged by scaling the amplitudes of the confidence values. 

Sixteen iterations (2* bit patterns) are performed to search for and invert the 
bits designated by the Golay decoder as the lowest four confidence bits in the 
codeword. The errors are counted for each iteration and the bit pattern with the 
Jowest cumulative confidence is selected as valid.7 

For both Golay decisions, there are 6 codewords (24 bits codeword x 6 
codewords superframe = 144 bits) generated for every superframe for the simu- 


latien, 


C. HAMMING CODES 

The Hamming (7.4) codes is a singie error correcting code where all combi- 
nations of two errors are detected: only one error may be both detected and cor- 
rected. This code was the first error detection and correction code discovered, in 


7 Bauds 10, 00. 01, and 11 correspond to phase changes 45, 135, -135, and -45 degrees, re- 
spectively. 
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Figure 5. Confidence Constellation [Ref. 2] 


1949 by R.W. Hamming (Ref. 18: p. 13-29}. The modified (8,4) Hamming code 
is used for this analysis, where, 1 parity bit is added to maintain code rate 
R=1/2. 

For the simulation, there are 18 Hamming codewords (8 bits/codeword x 18 


codewords = 144 bits/superframe) generated for every superframe. 


D. QUADRATIC RESIDUE CODE (QRC) 
Quadratic residues are used to specify the roots of the code generator 


polynomial [Ref. 19: pp. 92,93]. They are defined by the numbers, 


12. 22. 3? p-1\ 
> ’ | ee ee hes 2 


which are mod p reduced, where p is an odd prime number. Nonresidues are the 


numbers in the quadratic group that are not included by the residues. For ex- 
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ample, if p = 19, this would include all the mod 19 residues for 17, 27, 37,..., 


9°. The following table, Table 2 on page 17, shows these values. 


Table 2. QUADRATIC RESIDUES FOR 
P=19 


P| residues | non-resid 





The QRC investigated by GTE is the (48.24) code, with p = 47. This code 
is capable of detecting and correcting 100% of 5 (or fewer) bit errors and 62% 0 
of 6-bit errors. The probability of correcting 6-bit errors is calculated with com- 


binatorial logic as follows, 
(47.24) + 1 parity bit = (48,24) 
The number of 5 or less correctable errors is: 
C(4A7.5) + C(47-4) + C4723) + C(47,2) + C(47,1) = 1,729,647 


Also, the codeword with one error must be accounted for, C(47,0) = 1. Since 


there are 2-3 possible parity vector values (8.388.608), 


8.388.608 — 1.729.647 — 1 = 6.658.900 





is the number of possible parity comibinations remaining to map into 6 or more 
error combinations. Next, the number of 6-bit errors possible within the 48 bit 
codeword is: 


C(47,6) = 10,737,573. 
Therefore, 
6,658,960 + 10,737,573 = .6202 


This implies that the QRC will correct the sixth error in a 48 bit codeword 
62.02% of the time. 

There are 3 QRC codewords (48 bits,codeword x 3 codewords = 144 
bits'superframe) generated for every superframe in the simulation. 

The implementation of these EDAC schemes into a transmission simulation 
model will be discussed in Chapter IV. Also, the testing design approach and 
channel simulation models are introduced the next chapter. 





IV. TESTS AND SIMULATION 


A. DESCRIPTION OF THE PROGRAMS PROVIDED 
1. Introduction 
The program used in this analysis was provided by DIRNSA;R556, and 
Was previously used by DIRNSA to evaluate the use of Golay (hard and soft 
decision) coding for the protection of the CELP information bits. A separate 
subroutine of the Hamming code was also provided, for integration into the main 
program. Program subroutines neatly divided most of the program chores (i.e. 
bit generation, encoding, modulation, decoding, bit error count, etc.). 
2. Modifications 
The program code was loaded and run on a VAX 11/785 mainframe in 
the ECE Department at the Naval Postgraduate School after conversion from 
Sun I]] FORTRAN to Berkley FORTRAN 77. The program was further mod- 
ified to incorporate the Hamming code Simulation and a bit error counter for the 
QRC. The encode-decode functions of the QRC are not included in the simu- 
lation - only a counter based on the code’s probabilitv of error detection and 


correction within a codeword.8 


B. SIMULATION DESCRIPTION 

Figure 6 outlines the simulation process in block fashion. The process begins 
as a pseudorandom (PN) bit stream, representing the MI bits, is produced for 
each superframe. Next, the 72 MI bits are encoded by the chosen EDAC, with 
“no code” representing the default control. The encoder generates the 72 parity 
bits to formulate a completed codeword of 144 bits. 

The !44 bits are then hashed. This hashing process is in reality a pre- 
interleaving step that enhances the spread of the codeword bits throughout the 

8 At this wnting GTE 1s looking at simulating the QRC. However, the algorithmic imple- 


mentation is much more complicated than even the Golay codes. Therefore, actual QRC encoding 
and decoding is not attempted in this analysis. 








superframe.’ For each codeword the first bit is stacked, next the second, the third 
and so on. These realigned bits are then interleaved in Scheme II fashion, as 
described in Chapter II. To fill the rest of the 3360 bit superframe, an arbitrary 
PN bit stream is produced. The codeword bits are interleaved among these fill 
bits. 


SIMULATION BLOCK DIAGRAM 


72 bits - 


fillbits - 


*used for soft decision golay only 
**using either Golay (hard or soft decision), Hamming, QRC, or no code 





Figure 6. Simulation Block Diagram 
{ 


Before modulation and transmission, the bits are changed to baud using a 
dibit representation as shown in Table 3 on page 21. The entire superframe is 


modulated and transmitted, then sent over one of three user-chosen simulation 


9 This process is not reflected in the codeword interleaving table in Appendix C. 
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channels. The transmitted bits, now with simulated errors, are received and de- 
modulated, deinterleaved, stripped of the 3216 non-applicable bits, de-hashed, 
decoded and checked for errors. 


Table 3. BIT/BAUD EQUIVALENT 





Figure 7 illustrates the simulation design approach. Each section of the ma- 
trix represents one simulation run, for a total of 80 runs (5 codes x (1 
INMARSAT CH + 8 AWGN CH + 7 constant burst CH)). A simulation run 
consists of 200 consecutive superframes. In all cases, the first superframe was 
expended for synchronization purposes. Thus, there are 199 x 420 ms or 83.58 
seconds of established continuous information flow, which is likely much greater 
than an average message transmission. During the 83.58 second period, conti- 
nuity of the digital signal is maintained unless bit synchronization is lost due to 


fade or noise. 
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Figure 7. Simulation Design Approach 


All simulation runs record the number of errors (within the 144 EDAC 
bits/superframe - location of an error is unknown) out of the total number of 
EDAC bits processed. Also, the number of successful and failed (S/F) super- 
frames are recorded. Even one uncorrected error results in a superframe failure. 


C. CHANNEL SIMULATIONS 

‘the “mode” and “vary” variables provided for each channel are simulation 
selectable variables which determine the channel simulation mode. 

For the simulation model, during fading or burst mode sequence the signal 
was subjected to -24 dB S/N and the signal strength divided by 200. This was 
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injected to ensure a 50% Bit Error Rate (BER) during burst error as required bv 


the INMARSAT model. During non-burst error intervals, the S;N could be at 
any level. This parameter is user specified. For the two burst channel simu- 
lations, the S, N was set equal to 99.0 dB to ensure a clean signal. Therefore, a// 
of the errors are a result of the signal fade. In the Gaussian channel the S/N 
parameter was varied from 0 to 7 dB. For the Gaussian channel, simulation er- 
rors Were a result of poor S,N levels only. 
1. INMARSAT Burst (mode = I,vary = 0) 

For the random burst noise (signal fade) channel simulation, the “simpli- 
fied Land Mobile Radio channel model for IMARSAT Codes Evaluation Pro- 
posal (ICEP)” was used [Ref. 20: pp. 9,10]. The model is described as follows: 


“The model is derived from realistic empirical propagation measurements such that 
performance of voice codecs can be tested in the laboratory with replicable situ- 
ations. 90° is chosen as the link availability figure in the simplified model. The 
model is basically an ON-OFF model in which fading is in a binary state. The 
ON-state should correspond to no transmission whereas the OF F-state corresponds 
to 50°. BER.” 

The discrete probability density function (PDF) for the permissible noise 
burst widths are shown in Table 4 on page 24. The distribution of the burst 
widths for the INMARSAT channel is base on a random process. Therefore, the 
number of bursts for each burst width varied. For all but 200 msec burst widths, 
the actual number of bursts for each burst width were not significantly different 


from that predicted by the PDF. 
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Table 4. SIMPLIFIED 
FADE DU- 
RATION 
PDF 





For the burst-noise channel simulation, one 200 superframe run was per- 
formed for each of the four codes and one run with no coding. An example of 
3-second Golay-encoded signals, transmitted over an INMARSAT burst-noise 
channel, is shown in Figure 8 on page 25. 

2. AWGN Channel (mode = 0,vary = 0) 

For the additive white Gaussian noise (AWGN) channel simulation, 
signal-to-noise ratio (S/N) was varied from 0 to 7.0 decibels (dB) in 1.0 dB in- 
crements. Thus forty simulation runs (5 code modes x 8 S;N levels) were per- 
formed. 

3. Constant Burst Width (mode = I,vary = 1) 

For the constant burst width channel simulation, a constant fade depth 
(or burst Width) was inserted at random time intervals. Seven discrete fade 
widths -- 5, 10, 20, 40, 100, 200, and 400 milliseconds -- were run for each code. 
Thus, thirty-five simulation runs (5 codes x 7 burst widths) were performed, each 
having 200 superframes. The 5 and the 400 msec burst widths, not included in 
the INMARSAT model, were added to observe the limits of the model and the 


performance of the codes. 
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Figure 8. INMARSAT Channel Performance [Ref. 2] 


Figure 9 shows constant amplitude signals subjected to varying fade 
widths. One second translates to 2.38 superframes. 

Using the same outline as Chapter IV, the test results are presented in 
Chapter V by channel simulation type. The BER results and S/F resu!ts for each 


simulated code over each of the three channels is displayed and compared. 
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VARYING FADE WIDTH SIGNALS (1 SEC) 


burst width .ms 





Figure 9. Varying Fade Width Signals [Ref. 2] 
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V. COMPARISON AND ANALYSIS OF TESTS 


For each of the three simulated noise channels, the codes are rated in terms 
of BER performance and in terms of number or percentage of successful (or 
“error-free”) superframe transmissions. BER results are expressed as the BER 
exponent.!@ This number is calculated by taking the log,, of the number of errors 
counted divided by the number of bits transmitted. 

The code abbreviations used in data tables are: 

® g(s) - soft decision Golay (24.12) 

® g(h) - hard decision Golay (24,12) 

@ none - no channel coding 

¢ ham - Hamming (8.4) 

® QRC - quadratic residue code (48,24) 


A. INMARSAT BURST CHANNEL 
1. BER Results 
Table 5 on page 27, reflects BER results as the log of the quotient 
(number of errors counted divided by the number of bits transmitted). All BERs 
are greater than 10°! for this simulation, Which implies greater than 1| error out 
of 100 bits. These performance values are intuitively reasonable for LMR oper- 


ations in an urban environment. 


Table 5.) A FUNCTION OF INMARSAT BURST NOISE, ALL CODES 


[few [etn none [ham | QRC_ 





Also, the difference between the hard and soft Golay is only about 0.5°%, 
while the QRC cutperforms the soft Golay by more than 1.6%. This mav be a 


reflection of the number of “long” bursts. For the pseudo-random process, the 


10 The convention in [Ret. 2] was followed. 
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g(s) suffered eight 200 msec fades while the QRC was only foiced to process 2. 
If these additional bursts happened to span the critical MI bits in the g(s) 
codewords, the BER would have increased.!! Looking at the INMARSAT PDF, 
Table 4 page 24, the 200 msec bursts only have a 1% probability of occurrance. 
The difference between the number of 200 msec fades for the g(s) and the QRC 
runs is therefore, insignificant. 

Figure 10 on page 29, derived from table 5, shows the BER performance 
advantage that the QRC has over the other codes for this simulation. The QRC 
detected and corrected 1.5% more errors than the g(s). 

2. Run Success/Failure Results 

Table 6 on page 28 reports the number of successful runs in the simu- 
lated INMARSAT channel, for each code.!2 While BER performance favors the 
QRC in an INMARSAT nvise channel, the QRC only ran an additional two 
successful superframes out of 199. That is less than a 1% performance differ- 
ence. Nevertheless. the QRC and g(s) are noticeably superior than the other 


codes and the uncoded channel. 


Table 6. INMARSAT BURST - #SUCCESSES/#FAILURES 


Ti) [at [none [ham [ORC 


FSuccesses 


#failures 


#Htotal runs 


% successful 





11 This simulation did not allow the control of the number of bursts of each width. Therefore, 
it was impossible to inject the codes with the same error occurances. 


12 If the number of successes and the number of failures do not add up to 199, the runs not 
accounted for are out-of-svynchronization (OOS) superframes. If there are too many errors the 
transmission falls OOS and retries until a success or a failure 1s achieved. In any of the channel 
sumulations the OOS runs are not counted toward the total number runs. 
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Figure 10. INMARSAT BER Performance 








3. Additional Tests 

Further testing was performed to check the validity of the INMARSAT 
Simulation results. An additional 1000 superframes were run for both the g(s) 
and the QRC. Again, the QRC outperformed the g(s) in BER by .62°%. The 
QRC also provided 2% more successful superframe runs than the QRC. 

There seems to be a trend, in that the performance difference between the 
two codes appears to narrow with more superframe runs. This hypothesis could 
be tested with even more superframe runs. Nevertheless, the BER difference of 
less than 1% leads to the conclusion that there is no significant difference be- 


tween the INMARSAT channel performance of the two codes -- g(s) and QRC. 


B. GAUSSIAN NOISE CHANNEL 
1. BER Results 
Table 7 on page 31, tabulates BER exponent for the 40, 200 superframe 
runs performed in the simulated Gaussian noise channel.!3 Figure 11 on page 32 
is derived from Table 7. and shows that g(s) outperforms QRC at S'N levels 
greater than 5 dB. The figure also shows that the g(h) outperforms the QRC at 
S N levels greater than 6 dB. At S,N levels less than or equal to 5 dB, the gis) 


offers a slight improvement over that of the QRC and the g(h).!4 


13 No errors were counted at 7 dB for g{s) out of codeword bits; the BER exponent of -4.0 
was used as a graph limit. 


14 For the Hamming code, extra data points were taken; thev are reflected in the graph only 
and not the table. 
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Table 7. BER AS A FUNCTION OF GAUSSIAN NOISE LEVEL, ALL CODES 


-0.6055 -0.6421 -0.654 
-0.7447 -0.719 -0.7525 
-0.9547 -0.8153 -0.9014 


-1.2388 -0.9281 -1.1537 
-1.6126 -1.0605 -1.556 

~2.2676 -1.2269 -2.1586 
-3.2343 -1.4056 -2.6198 
-4.0 - 1.6234 -2.8339 





Figure 12 on page 33 is another representation of the same data. Code 
performance is compared at each S,N level measured. For each code, as the S'N 
increases the BER decreases. There are no surprising trends for any code. These 
two factors -- S N and BER -- are inversely proportional and the comparison 
chart best illustrates this response. Also, the difference between the codes re- 
sponse at Various S.N levels is more easily observed on this comparison chart 
than on the previous line graph. 

2. Run Success/Failure Results 

Table 8 on page 34 displays the number and the percent of successful 
runs with the OOS runs excluded from the total. These calculations were made 
using the total number of successful runs over the total number of runs per- 
formed. The g(s) code outperforms other codes in the gaussian environment, es- 


pecially at lower S N, which are more realistic in an LMR environment. 
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Figure 12. Gaussian BER Performance (bar chart) 








Table 8. #SUCCESSES/%SUCCESSFUL RUNS VS. GAUSSIAN NOISE 
LEVEL, ALL CODES 


| dB | es) #% | eth) #% | none 4% | ham #i% 


0.0 0;0 0;0 0/0 
2:10 0/0 0/0 0:0 
19 9.6 3/1.5 i 0;0 7;3.6 

62, 31.5 26/13.2 2/1.3 35:17.6 


122;61.3 72:36.9 7/3.6 107/53.8 
175,87.9 140.71.0 41/20.6 177;88.9 
194.97.5 186,93.5 92/46.2 19497.5 
197,99.0 194.97.5 149 74.9 198,99.5 





The line graph, Figure 13 on page 35 compares the percent success of 
each code in the Gaussain noise channel. For this fiqure of merit, the g(s) code 
outperforms both the g(h) and QRC, especially at S’N levels between 2-5 dB. 
The QRC outperforms g(h) at S-N levels between 2-7 dB. 

The stacked bar comparison graph, Figure 14 on page 36 compares the 
overall success of each code, combining all SN level runs. For total number of 


runs, g(s) out performs QRC by more than 50 successes. 


C. CONSTANT BURST WIDTH CHANNEL 
1. BER Results 
Table 9 on page 37 represents BER exponents for 35, 200 superframe 
runs. Figure 15 on page 38 shows the BER performance recorded in table 9, 
with the area of critical code performance highlighted on the graph.!5 The g(s) 


13 from Table 4 page 21), INMARSAT Noise Burst Width: P(10 msec burst) = 0.80; 
P( 10-40 msec burst) = 0.95. 
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code outperforms the QRC at burst widths from 10-40 msec; the advantage di- 


minishes at burst widths greater thatn 40 msecs. 


Table 9. BER AS A FUNCTION OF CONSTANT BURST WIDTHS, ALL 
CODES 


“4, 
-3.1759 
-2.0301 


-1.3516 
-0.7595 
-0.585 

-0.4789 





Figure 16 on page 39 also displays the effect that the varving burst 
widths have on each code. At the 20 msec burst width all three of the codes -- 
g(s), QRC, and g(h) -- begin to display a similar performance. For a 100 msec 
burst width condition, a// codes perform equally poorly. Code response trends -- 
a consistant increase in the BER exponent as the burst width increases -- are as 
expected. 

2. Run Success/Failure Results 

Table 10 on page 40 displays the number and percent successful runs as 
function of constant burst widths, for all codes tested. At narrower levels of burst 
width, the g(s) code outperforms the other coded or uncoded channels. The ad- 
vantage becomes less significant as the width of noise burst excecds 20 msecs. 


Again, the g(s) has the overall advantage of 4.9% more successful runs. 
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Figure 15, 
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Table 10. #SUCCESSES/%SUCCESSFUL RUNS VS CONSTANT BURST 
WIDTHS, ALL CODES 


widths 


197;99.5 196:99.5 98/49.2 187;94.0 198,100 
197;99.0 188;94.5 93/46.7 152/76.4 195;98.0 
178;89.4 158/79.4 85/42.7 118/59.3 166, 83.4 
139.69.8 97/48.7 74'37.6 77/38.7 107,/53.8 


6633.2 4824.1 62;31.5 48 24.1 41/20.6 
39, 19.6 19/425 39;20.0 26;13.1 2613.1 
1/5.9 0,0 8/30.0 28, 14.1 19 11.0 


eo 


The line graph, Figure 17 on page 41, compares the relative number of 





successes for each code. A more robust code with greater interleaver depth may 
be needed to satisfy performance objectives at these burst widths.!6 The g(s) code 
performs best, but it’s percent success degrades as burst widths exceed 40 msecs. 
It is interesting to note that the uncoded channel surpasses all codes in perform- 
ance except the g(s) for burst widths of 100 and 200 msecs. 

Figure 1&8 on page 42 also compares the percent success for each code. 
Again, the performance advantage of the g(s) and QR codes degrades beyond a 
40 msec burst to the 100 msec burst. Figure 19 on page 43 shows the perform- 
ance trends as burst widths lengthen. With the exception of the g(s) and QRC, 
all codes achieve less than 50° success at burst widths greater than or equal to 


40 msec. 


16 More analysis would needed to reach a decisive conclusion. 
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Figure 17. 
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Figure 18. 


42 
































E 
+2) 
= 
3 
3 
zg 
3 
rs) 
= 
mee fae Gee NS eR ey op VLLLLE 
en : ERR RRR RR ee ODP a 
: ‘ ita! atatatasat ee ty © 
: Ciba pis Littl. 
ee ae a a a Oe ee eS Ce 
BAO Migs ae raat cs ae ee ae Seca We SE ieee 
a M erateeeie'e OT E 
8 £ : ‘ : : : . ‘GGutt eicetee ei 
no: ODI EI 
” bad . cee Lette tasers meson tee teem crete ererensecestectereseesesecsccsecreenesaersscorsetmcetem ated 
ig See eee eg a, ee 
Oo § Rie renee was Sa Pan ee Tat ae? MCT ee Ve Tae See ee edn Sy. 
om ee OR ON ae es —bhttitita 
a ITE 0-0-0 0-0-00 0.6 -0-0:0°00.0.0.00 aS 
eae . PRE Sg 
ee ee a LLLILEL LAN ANA, 
° 
° 


SSSI 


: : warererere wererevavaverer “sO eLrOTAAOnOCATOTOLONS. = 
fsssscomasuinnessssuiuuuia unui HEH 


VLILLLLLLL LLL ALAA LALLA ELLE 2222 TISTEMA: 





: TREE 


als) 





. . . x? 
- RED seem zi Soe rree perro 
: LILLIE ULLEDUILIDIDOLEDILDPIDILITILIDIEA 








Figure 19. Constant Burst Width - S/F (comparison) 
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D. COMPARISON OF ERROR POSSIBILITIES 

Given the larger codeword, by probablistic calculations, the QRC leads in 
error correction capability since there are a greater number of possible errors 
corrected. See Figure 20 on page 45 for an understanding of the possibilty of 6 
errors corrected out of 48 bits for each code. Notice that the Golay calculation 
considers the hard decision number of possible errors only. The soft-decision 
Golav code yielded a 5% better success rate for the INMARSAT channel (see 
Table 6 page 28). It was also observed that the number of possible errors cor- 
rected is much closer to the QRC calcution - up to seven errors can be corrected 
per codeword. 

In Chapter VI, the preceeding results will be summarized and recommen- 


dations presented. 








COMBINATIONS of 6 errors/48 bits 
corrected 


62% x P (C(48,6),1) = 7,610,512 


a ee 


P(C(24,3),2) + C(24,3)C(2,2,)P(C(24,3),0) = 4,096,576 


Golay | | | 


24 bits 


P(C(8,1),6) + C(8,1)C(6,2)P(7,4) + ... + C(8,1)C(6,6)P(7,0) =159,944 


Hamming 


8 bits 





Figure 20. Possible Errors Corrected 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER 
STUDY 


A. CONCLUSIONS 

Three simulated noise channels -- INMARSAT random burst, AWGN, and 
constant burst -- were used to analyze four EDAC codes; Golay hard decision, 
Golay soft decision, Hamming and QRC. Before the study commenced, the 
INMARSAT random burst channel model was determined to be the most real- 
istic type of LMR noise channel. It was believed by NSA that the soft Golay and 
QRC codes would be very close perfomance competitors. The results obtained 
show that the soft decision Golay slightly outperforms the QRC based on the 
following: 


e For the Gaussian channel, the BER exponent for the g(s) is significantly less 
than that of the QRC at S'N of 5 dB and below. The g(s) code also had 
greater percentage of successful runs than the QRC for S,N levels from 2-4 
dB. 


e For the Constant Burst Width channel, the g(s) marginally outperforms the 
QRC when subjected to 10 msec burst noise. For each code, percentage of 
successful runs is similar at this burst width; the g(s) achieves 1%o greater 
success. Over all burst widths, the g(s) achieves 4.9% more successful 
superframe runs. 


e For the INMARSAT random burst channel the results are very close. The 
QRC outperformed the g(s) by detecting and correcting 1.5% more errors 
and running almost 1% more successful superframes. For the additional 
1000 superframe runs, the results are even closer. Since the performance 
margin is so close, the conclusion is that there is no significant difference 
between the two codes for this simulation. More superframe runs may fur- 
ther support this conclusion. 


The margin of performance is very close, as expected, between the soft deci- 
sion Golay and the QRC. The biggest disadvantage in QRC implementation 
versus soft decision Golay is the difficulty of algorithmic implementation and the 
larger number of codeword and parity combinations to chose from. This trans- 


lates into greater processing time and greater memory storage. 








Based on these two considerations, soft decision Golay is recommended as the 
best code for FS 1024. 


B. RECOMMENDATIONS FOR FURTHER STUDY 
1. Study Continuation 

From experience it was learned that a more concentrated level of effort 
should have been paid to the 10-40 msec burst width region (95% of the burst 
fades for the INMARSAT PDF). More superframes should be subjected to these 
burst widths to provide more conclusive findings. 

Also, additional INMARSAT channel iterations along with a statistical 
analvsis!7 may shed light on the difference that the long burst widths may have 
on the performance of the g(s) versus the QRC. 

2. Implementation 

Another consideration and topic for further study is the feasibility and 
architecture of hardware implementation of candidate EDAC algorithms. Se- 
veral Digital Signal Proccessor (DSP) Integrated Circuits (ICs), such as the 
DSP32C and the TMS320C30, are being considered for hardware implementa- 
tion. Microprocessor (assembly language) code of the soft decision Golay is al- 
ready under development » DIRNSA.18 Since the QRC codebook is much larger 
than the soft decision Golay, IC memory (itself limited by chip size) may be a 
limiting factor for hardware implementation of the QRC. 

3. Other Codes 

The scope of the thesis was narrowed to the four block codes: however, 
Other possibilities should be considered in future studies. Simulation and testing 
of convolutional, other linear block, and concatenated codes may result in the 
discovery of an even better code for LMR application. A comparison of appro- 
priate tradeoffs -- such as processing time, difficulty of algorithmic implementa- 


17 The program must be modified to obtain the necessary statistical information. such as - the 
bit position of the error, position of the fade, the number of errors each fade caused, etc. 


18 This EDAC was chosen as the code to protect the digitized voice (CELP processed) infor- 
mation bits. 
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tion, and code performance -- should be weighed in order to fairly judge the best 
code for the MI bits. 
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APPENDIX A: FORTRAN CODE 


A. MAIN DRIVER 
c 
c 
c 
¢ 


COCO CO OIG Gm iol OR mii iii i ai gk a dokok ak kk 


14*240=3360 bits per superframe (sf) 
72+72=144 bits per mi 
3360-144=3216 bits of fill per sf 
3360/2=1680 bauds per sf 


This program generates a superframme of 3360 bits. 

72 bits are of interest in the detection and correction 

of errors. Three different EDAC schemes simmulated here 
for comparison of the best scheme under various conditions. 
The three codes are Golay (24,12) using soft and hard 
decision logic, Hamming (8,4), and Quadradic Residue Code 
(QRC) (48,24). Also, a run with no coding is performed as 
a control case. 


The entire superframe is interleaved, transmitted over a 
simulated AWGN channel, received, deinterleaved, and 
finally - checked for errors. The run in considered a 
failure if there are ANY errors. 


Options for the simmulated channel include variance of the 
s/n during non-fade (the fade s/n was held at -24.0 dB), and 
two burst modes. The first burst mode used the INMARSAT 
values of varying length from 10-200ms for LMR. The second 
burst mode allows for a constant burst length input. 


CORO CK OOO OR GO Ok kok ik kok iok kk ak 


program mainedac 


CAO eR dD @ C1 a7 at 1 O08 IO a OR IOI GIO 


integer*4 tbaud(1680) 

integer*4 fillbits(3216) 

integer*4 parity(72) 

integer*4 rbaud(1680) 

integer*4 rbcnt,tbcnt ,decodcnt1,decodcnt2 
integer*4 modegolay,ham,parnum,qrc,sf,nummi 
integer*4 index,prevoutdibit ,table(16) 
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integer*4 tbits,bits,txbits,rbits 

real xdisp,ydisp,confh,confl 

real samples(3) ,noise,noisef ,bprob,probq 

integer*4 mode,iseed,vary 

integer*4 testcnt,tbit (3360) ,rbit(3360) ,intlvtable(144) 
integer*4 mibits(72) ,codeword(8) ,hmatrix(8) ,syndrometable(8) 
integer*4 testcnt ,paritybit 

real intlvconf (3360) ,rcodeconf (144) ,rhashconf (144) 

real rfillconf (3216) 

real conf (24) ,gaus(256) 

integer*4 code(144) ,hashcode(144) ,qrcbits(144) ,hambits (72) 
integer*4 hashtable(144) 

integer*4 rcode(144) 

integer*4 itoterr,qrccwerr,qrcsferr,bwidth ,myerror,hamerr 
integer*4 iallerr, itotbitct 

logical insync 

integer*4 success,bits,fail,nerror,intlvrcnt,ifillrcnt,ihashrcnt 
integer*4 rhashcode(144) ,rfillbits (3216) 

integer*4 isym(2047) 

integer*4 idata(24) ,kdata(24) ,ierr(5) ,qdata(48) ,hdata(4) 
common /b1kS/idata 

common /b1k7/conf,isym,kdata,ierr 


Ca OOO A OOaREDATA and TABLES Sado Goa oon a Ia oR Aka 


et ce a oe oS a oS OO 





myerror = 0 

data init,initfill /1,1/ 

data table /2,0,3,1,3,2,1,0,0,1,2,3,1,3,0,2/ 

if ((ham.eq.0).and.(qrc.eq.0)) then 

data hashtable / 

1,25,49,73,97,121,2,26,50,74,98,122, 
3,27,51,75,99,123,4,28,52,76,100,124, 
§,29,53,77,101,125,6,30,54,78,102,126, 
7,31,55,79,103 ,127,8,32,56,80,104,128, 
9,33,57,81,105,129,10,34,58,82,106,130, 
11,35,59,83,107,131,12,36,60,84,108,132, 
13,37 ,60,85,109,133, 
14,38 ,61,86,110,134, 
15,39,62,87,111,135, 
16 ,40,63,88,112,136, 
17,41,64,89,113,137, 
18 ,42,65,90,114,138, 
19,43,66,91,115,139, ° 
20 ,44,67,92,116,140, 
21,45,68,93,117,141, 
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20 


25 


eS a 


on 





22,46 ,69,94,118,142, 

23 ,47,70,95,119,143, 

24,48,71,96,120,144/ 
end if 


if (qrc.eq.i) then 
data hashtable / 

1,49,97,2,50,98,3,51,99,4,52,100, 
§,53,101,6,54,102,7,55,103,8,56,104, 
9,57,105,10,58,106,11,59,107,12,60,108, 
13,61,109,14,62,110,15,63,111,16,64,112, 
17,65,113,18,66,114,19,67,115,20,68,116, 
21,69,117,22,70,118,23,71,119,24,72,120, 
25,73,121,26,74,122,27,75,123,28,76,124, 
29,77,125,30,78,126,31,79,127,32,80,128, 
33 ,81,129,34,82,130,35,83,131,36,84, 132, 
37 ,85,133,38,86,134,39,87,135,40,88,136, 
41,89,137 ,42,90,138,43,91,139,44,92,140, 
45 ,93,141,46,94,142,47,95,143,48,96,144/ 

end if 


if (ham.eq.1) then 
data hashtable / 

1,9,:.7,25,33,41,49,57,65,73,81,89,97,105,113,121,129,137, 
2,10,18,26,34,42,50,58,66,74,82,90,98,106,114,122,130,138, 
3,11,19,27,35,43,51,59,67,75,83,91,99,107,115,123,131,139, 
4,12,20,28,36,44,52,60,68,76,84,92,100,108,116,124,132,140, 
§,13,21,29,37,45,53,61,69,77,85,93,101,109,117,125,133,141, 
6,14,22,30,38,46 ,54,62,70,78,86,94,102,110,118,126, 134,142, 
7,15,23,31,29,47 ,55,63,71,79,87,95,103,111,119,127,135,143, 
8,16,24,32,30,47,56 ,64,72,80,88,96,104,112,120,128,136,144/ 

end if 


read in table for golay decoder from disk file 

open (3,file=’isytab.dat’ ,status=’old’ ,form=’ formatted’) 
read (3,20)isym 

format (10(16) ) 

close(3) 

read in gaussian distribution from disk file 

open (4,file=’gaussian.dat’ ,status=’old’ ,form=’formatted’ ) 


read (4,25)gaus 
format (8(£6.0,1x)) 
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close(4) 


CaaeeeE Tn teractive Tnputaeiiesiieaic a OO OO IOC 


write(6,*)’how many 420ms superframes would you like to run?’ 
read(S,*) sf 
write(6,*)’what golay mode ? (1=soft,2=hard,3=nogolay) = 7?’ 
read (5,*) modegolay 
if (modegolay.eq.3)then 
write(6,*)’would you like to try the Hamming code? 
*(1*yes,0=no) .. .’ 
read(5,*) ham 
if(ham.eq.1) qrc=0 
if (ham.eq.0) then 
write(6,*)’would you like to perform a QRC error count? 
*(1=yes,0=no) ... ’ 
read(5,*) qrc 


end if 
else 
ham=0 
qrc=0 
end if 
write(6,*) ’what is s/n ratio in dB during fade = ? (real #)’ 
c should be -24dB 
read (5,*) noisef 
write(6,*) ’what is s/n ratio in dB during no fade = ? (real #)’ : 


read (5,*) noise 
write(6,*) ’what is mode ? (0 for non-fading, 1 for fading) 
read (5,*) mode 
if (mode.eq.1)then 
write(6,*)’what is the time seed? (any positive integer) 
read (5,*) iseed 
write(6,*)’would you like to vary burst prob and width? 
*(izyes,0=no) ...’ 
read (5,*) vary 
if(vary.eq.1) then 
write(6,*)’what burst prob?(r .94056 INMARSAT) = ’ 
read (5,*) bprob 
write(6,*)’what burst width?(i 120=10ms, INMARSAT) = ° 
read (5,*) bwidth 
end if 
else 
iseed = 1 : 
end if 


it} 
- 
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cxeeee% write the superframe interleave table for the 144 codeword bits 
ceeeeee only throughout frames 4 through 14 (frame 4 starts at bit 721) 


do 35 m=1,9 
do 30 mm=i1,8 
intlvtable(((m-1)#16)+((mm-1)*2)+1)= 
1 721+ ((m-1)*240)+( (mm-1)*30) 
intlvtable(((m-1) *16)+((mm-1) *2)+2)= 
1 721+( (m-1)*240)+( ( (mm-1) *30) +10) 
30 continue 
35 continue 


Cee Tri tia lize GIO GOS IA OR #08 


itoterr = 0 
qrcsferr = 0 
qrcecwent = 0 
success = 0 


fail = 0 
ifirst=1 
tbcnt=1680 


rbcnt=1680-10 


cee*#* Begin for sf# of sUpeT frames ee H HAE OOO IA 
c#**ee* don’t start counting results until n=about 3 so that error counters 
c*¥**** are in sync (test if insync and insyncfill .eq..true) 


do 500 n=i,sf 
c loop for 1680 baud/sf - generate, tx, sim, rcv, etc. 
do 400 mk=1,1680 
tbcnt=tbent+1 
if (tbcnt.eq.1681) then 
tbent=1 


ceeeke FIL] mibits 
if (qrc.eq.1)nummi=72 
if (qrc.eq.0)nummi=36 
do 40 m=1,nummi 
call bitgen11(2,init,tbits) 
init=0 
if ((qrce.eq.0).and. (ham.eq.0)) then 
mibits(((m~1)*2)+1)=0 
mibits(((m-1)*2)+2)=0 
if ((tbits.eq.1).or. (tbits.eq.3))then 
mibits(((m-1)*2)+1)=1 
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end if 


if ((tbits.eq.2).or.(tbits.eq.3))then 
mibits(((m-1)*2)+2)=1 
end if 
end if 
if(ham.eq.1) then 
hambits(((m-1)*2)+1)=0 
hambits(((m-1)#2)+2)=0 
if ((tbits.eq.1).or.(tbits.eq.3))then 
hambits(((m-1)*2)+1)=1 
end if 
if ((tbits.eq.2).or. (tbits.eq.3))then 
hambits(((m-1)*2)+2)=1 
end if 
end if 
if(qrc.eq.1) then 
qrcbits(((m-1)*2)+1)=0 
qrcbits(((m-1)*2)+2)=0 
if((tbits.eq.1).or. (tbits.eq.3))then 
qrcbits(((m-1)*2)+1)=1 
end if 
if ((tbits.eq.2).or.(tbits.eq.3))then 
qrcbits(((m-1)*2)+2)=1 
end if 
end if 
40 continue 
if (qrce.eq.1) go to 95 


ce*eeee fil] parbits 
if (ham.eq.1) parnum=18 
if (Cham.eq.0).and.(qrc.eq.0)) parnum=6 
call matrixgen(8,4,hmatrix,syndrometable) 


do 70 m=1,parnum 
cx**e* Dit encoder for hamming 
if (ham.eq.1)then 
do 48 mm=1,2 
codeword ( ((mm-1)*2)+1)=hambits(((m-1)*4) 
1 +(((mm~-1)*2)+1)) 
codeword ( ((mm-1)*2)+2)=hambits(((m-1) #4) 
1 +(((mm-1)*2)+2)) 
codeword ( ((mm-1)*2)+5)=0 
c codeword ( ((mm-1)*2)+6)=0 : 
48 continue 
call encodeham(8,4,hmatrix,paritybit , codeword) 
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c#¥**** load hamming codewords (parity + data) 30K 
do 50 mm=1,8 
code(((m-1) #8) +mm)=codeword (mm) 
50 continue 
else 


ce*xekk generate parity bauds for golay here 
if (qrc.eq.0)then 


do 55 mm=1,6 
idata(((mm-1)#*2)+1)=mibits(((m-1)*12) 
1 +(((mm-1)*2)+1)) 
idata(((mm-1)*2)+2)=mibits(((m-1)#*12) 
1 +(((mm-1)*2)+2)) 
55 continue 
end if 


call golenc 


do 60 mm=1,12 
parity (((m-1)*12)+mm)=idata(12+mm) 


60 continue 
end if 
70 continue 


cee*** load Golay and QRC codewords (bits) 2ERKH HE OEE Rito 
if ((ham.eq.0).and.(qrc.eq.0)) then 
do 90 m=1,6 
do 80 mm=1,12 
code(((m-1)*24)+mm)=mibits(((m-1)*12)+mm) 
code(((m-1)*24)+mm+12)=parity(((m-1)*12)+mm) 


80 continue 

90 continue 
end if 

95 continue 


if(qrc.eq.1) then 
do 97 m=1,144 
code(m)=qrcbits(m) 
97 continue 
end if 


c**x*exe* Scramble the codewords 
do 100 m=1,144 
hashcode(m)=code(hashtable(m) ) 
100 continue 
do 110 m=1,1608 


wii 
ur 














c get the rest of the superframe’s bits 
call bitgen11fill1(2,initfill,tbits) 
c same as bitgenii but a separate routine 


initfill=0 

fillbits(((m-1)*2)+1)=0 

fillbits(((m-1)*2)+2)=0 

if((tbits.eq.1).or. (tbits.eq.3))then 
fillbits(((m-1)*2)+1)=1 

end if 

if ((tbits.eq.2).or. (tbhits.eq.3))then 
fillbits(((m-1)*2)+2)=1 

end if 

110 continue 


c***** load the superframe with codewords and fill 
intlvcnt=1 
ihashcnt=1 
ifillcnt=1 
do 130 m=1,3360 
if (m.eq. (intlvtable(intlvcnt)))then 
tbit (m)=hashcode(ihashcnt) 
ihashcnt=ihashcnt+1 
intlvcnt=intlvcnt+1 
else 
tbit(m)=fillbits(ifillcnt) 
ifillent=ifillent+1 
end if 
130 continue 


c***k* bits to bauds 
do 140 m=1,1680 
if (tbit(((m-1)*2)+1) .eq.0)then 
if (tbit(((m-1)*2)+2) .eq.0)tbaud(m)=0 
if (tbit(((m-1)*2)+2) .eq.1)tbaud(m)=2 
else 
if (tbit(((m-1)*2)+2) .eq.0)tbaud(m)=1 
if (tbit(((m-1)*2)+2) .eq.1)tbaud(m)=3 
end if 
140 continue 
end if 


txbits=tbaud (mk) . 


c***ee* differential phase encoding 








index=(4*txbits)+prevoutdibit 
txbit3=table(index+1) 
prevoutdibit=txbits 


c¥*x*** transmit over simulated channel and receive 
call tx(txbits,samples) 
if(vary.eq.1)then 
call chsim(samples ,noise,noisef ,mode,iseed,gaus, 
1 bprob,bwidth) 
else 
call sim(samples ,noise,noisef ,mode,iseed,gaus ,n,sf,mk) 
end if 
call rcv(samples,rbits,xdisp,ydisp,confh,confl) 


c*x*e** The confidence calculations are used for softgolay only ***** 
rbcnt=rbcnt+1 
rbaud(rbcnt-1)=rbits 
intlvconf (((rbcnt-2)*2)+1)=confl 
intivconf (((rbcnt-2) *2)+2)=confh 

if (rbcnt.eq.1681) then 
testcnt=testcntt+i 
rbcnt=1 
if (ifirst.eq.1) then 
ifirst=0 
else 








cxk*ee* unpack and decode the bits and count any errors in the mi 
c this would be 18 calls to hamming or 6 calls to golay 
c if any of the 72 mi bits are in error then we fail and 
c*¥**** Dauds to bits 
do 150 m=1,1680 

rbit(((m-1)*2)+1)=0 

rbit(((m-1)*2)+2)=0 

if ((rbaud(m) .eq.1).or. (rbaud(m) .eq.3)) 


1 rbit(((m-1)#*2)+1)=1 
if ((rbaud(m) .eq.2).or. (rbaud(m) .eq.3)) 
1 rbit(((m-1)*2)+2)=1 
150 continue 


c*¥**e* deinterleave bits from the superframe 
intlvrent=1i 
ihashrent=1 
ifillrent=1 
do 160 m=1,3360 
if(m.eq. (intlvtable(intlvrcnt)))then 











rhashcode(ihashrent)=rbit (m) 
rhashconf (ihashrcnt)=intlvconf (m) 
ihashrcnt=ihashrcnt+1 
intlvrent=intlvrcnt+1 

else 
rfillbits(ifillrent)=rbit(m) : 
rfillconf (ifillrcnt)=intlvconf (m) 
ifillrcnt=ifillrent+i 

end if 

160 continue 


c****e* descramble the codeword bits 
do 170 m=1,144 
rcode(hashtable(m) )=rhashcode(m) 
rcodeconf (hashtable(m) )=rhashconf (m) 
170 continue 


c*#*ke decode the codewords (baud) to OOOO 
if (qrc.eq.1)then 
decodcnt1=3 
decodcnt2=24 
end if 
if (ham.eq.1)then 
decodcnti=18 
decodcnt2=2 
end if 
if ((qrc.eq.0).and. (ham.eq.0))then 
decodcnt1=6 
decodcnt2=6 
end if 
do 200 m=1,decodcnti 
if ((ham.eq.0).and. (qrc.eq.0)) then 
do 180 mm=1,24 
idata(mm)=rcode(((m-1)*24)+mm) 
conf (mm) =rcodeconf (((m-1)*24) +mm) 
180 continue 
end if 
if (ham.eq.i)then 
do 181 mm=1,8 
codeword (mm)=rcode(((m-1)*8) +mm) 
181 continue 


end if 
if (qrc.eq.1)then 
do 182 mm=1,48 
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qdata(mm)=rcode(((m-1)*48)+mm) 
182 continue 
end if 


if (modegolay.eq.1)call soft 
if (modegolay.eq.2)call hard 
if (ham.eq.1)then 
call decodeham(8,4,paritybit ,myerror ,hmatrix, 
& syndrometable,codeword) 


c¥*ee* strip off the first 4 bits of each codeword MI bits. *«*«+«#x*#«( 
do 183 mm=1,4 
hdata (mm) =codeword (mm) 
183 continue 


if (myerror.eq.i)then 
hamerr = hamerr + 2 
write(6,*)’decodeham detected errors ’ 
end if 
end if 





cx**** Count errors in 72 MI bits (golay=12/24,hamming=4/8) 
c or all 144 bits (qrc=48/48) 
cx**** before the errors are counted they are changed to baud 
CR OOOO GIOIA OOO OIG IGG RIG III tte 
qrccwerr=0 
do 190 k=1,decodcnt2 
if ((qre.eq.0) .and. (ham.eq.0))then 
if (idata(((k-1)*2)+1) .eq.0)then 
if (idata(((k-1)*2)+2) .eq.0)bits=0 
if (idata(((k-1)*2)+2) .eq.1)bits=2 
else 
if (idata(((k-1)*2)+2).eq.0)bits=1 
if (idata(((k-1)*2)+2) .eq.1)bits=3 
end if 
end if 
if (ham.eq.1)then 
if (hdata(((k-1)*2)+1).eq.0)then 
if (hdata(((k-1)*2)+2) .eq.0)bits=0 
if (hdata(((k-1)*2)+2) .eq.1)bits=2 
else 
if (ndata(((k-1)*2)+2) .eq.0)bits=1 
if (hdata(((k-1)*2)+2).eq.1)bits=3 
end if 


End if 





if (qre.eq.1)then 
if (qdata(((k-1)*2)+1) .eq.0)then 
if (qdata(((k-1)*2)+2) .eq.0)bits=0 
if (qdata(((k-1)*2)+2) .eq.1)bits=2 
else 
if (qdata(((k-1)*2)+2) .eq.0)bits=1 
if (qdata(((k-1)*2)+2) .eq.1)bits=3 
end if 
end if 
call biterriia(bits,2,insync,nerror) 
itoterr=itoterr+nerror 
if (qrc.eq.1)qrccwerr=qrccwerr+nerror 
190 continue 


c**ee* This is the QRC BER counter section 33GB OO IO OR Hoe 
c errors/codeword are checked (for decodcnti=3 cws) 
c QRC(48,24) will correct 5 errors/cw and 62.02% of the 6th errors 
COR OR RR Oi Gm ok OK OR ok koko kok kok gokotok dak ak koto ok kk 
if(qrc.eq.1) then 
if (qrccwerr.eq.0)then 
write(6,*)’there were no errors in QRC cw’ 
end if 
if ((qrcewerr.le.5).and. (qrccwerr.gt.0)) then 
write(6,*)’# errors corrected = ’,qrccwerr 
qrccwerr=0 
end if 7 
if (qrccwerr.eq.6)then 
probq = rand(0) 
if (probq.le. .6202)then 
write(6,*)’qre 6th error corrected’ 
qrccwerr=0 


end if 
end if 
qrcsferr=qrccwerr+qrcsferr 
end if 
195 continue 
200 continue 


c***** print error results 
c if any errors in total error count , the run is a failure 
CRISES IOI OIG IOI UIOI IO II UIOI IOI IOI OI IOI I IO GG I Kai 
if(qre.eq.i) then 
itoterr = qrcsferr 
iallerr = iallerr + qrcsferr 
itotbitct = itotbitct + 144 
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end if 
if (itoterr.eq.0)then 
success=success+1 
write(6,*)’success = ’,success 
end if 
if (itoterr.ne.0)then 
fail=fail+i 
write(6,*)’?fail = ’,fail 
write(6,*)’total errors in sf = ’,itoterr 
end if 


c*¥**** reinitialize for new superframe SiO OR ok kkk kok 

qrcesferr=0 

itoterr=0 

hamerr=0 

if(.not.insync)then 
fail=0 
success=0 
testcnt=0 
itotbitct=0 
iallerr=0 

end if 


c*#*ek* OF use the fill bits to measure ber in non-fading noise 


c call biterriiafill(bit,2,insyncfill,nerror) 
end if 
end if 
400 continue 


write(6,*) ’ you just completed superframe # ’,n 
if ((qrce.eq.1).and. (mod(n,50) .eq.0) )write(6,*) 


1 *jallerr= ’,iallerr,’ itotbitct= ’ ,itotbitct, 
1 n= ’,n 
500 continue 


c***e#e* print final counts and results 
write(6,*)’ the total number of failures = ’,fail 
write(6,*)’ the total number of successes = ’,success 
if (modegolay.eq.1) write(6,*)’ This was a soft golay run’ 
if (modegolay.eq.2) write(6,*)’ This was a hard golay run’ 
if (modegolay.eq.3) write(6,*)’ This was a nogolay run’ 
if (ham.eq.1)write(6,*)’ This was a Hamming run’ 
if (qrc.eq.1)write(6,*)’ This was a QRC BE ount run’ 
write(6,*)’ This run used fading (1=yes,0=no) ”°,mode 
if(vary.eq.i)write(6,*)’fading width, bwidth = ’,bwidth 


61 











write(6,*)’ s/n during fade (dB) = ’ ,noisef 
write(6,*)’ s/n during non-fade (dB) = ’,noise 


a eS Se OS SS SS OP ee er Se rr 


end 










B. SUBROUTINES 





TTTTItitirririritititirrirtrirriri titi rire rir irr rire etre res 
subroutine biterriia(data,nbits,insync,nerror) 


integer*4 data, nbits, nerror, rxdata, shindex, error 
integer*4 index, errcnt, bitcnt, i, itemp 
integer*4 nblck, btest, lshift, rshift, xor 
logical insync 
data errcnt, bitcnt / 0, 0 / 
data rxdata, index / 0, 1 / 
data nblck /100/ 


rxdata = rshift(rxdata,nbits) 
call mvbits(data,0,nbits,rxdata,11-nbits) 


shindex = lshift(index,2) 

itemp = xor(index,shindex) 

index = rshift(index,nbits) 

call mvbits(itemp,2,nbits, index, 11-nbits) 


nerror = 0 
error = xor(index,rxdata) 
if (error .ne. 0) then 
do 9 i =11-nbits, 10 
if (btest(error,i).eq.1) nerror = nerror +1 
9 continue 
- end if 


errcnt = errcnt + nerror 
bitcnt = bitcnt + nbits 


if (insync) then 
if (bitcnt .ge. nbJ-k) then 
if (float(errcent) .gt. .35 * float(bitcnt)) then 


insync = .false. 
nblck=100 
else 
write(6,*) ’errors in ’,nblck, ’ bits =’, errcnt 
nblck=3000 
c to span fades 
‘ errcnt = 0 
bitcnt = 0 
end if 
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end if 
end if 


if (.not. insync) then 
if (bitcnt .ge. 20) then 
if (float(errent) .le. .2 * float(bitcnt)) then 
insync = .true. 
write(6,*) ’insync = true’ 
else 
if (bitcnt.eq.20) write(6,*) ’not yet in sync’ 
index = rxdata 
if (index .eq. 0) index = 1 


end if 
errcnt = 0 
bitcnt = 0 
end if 
end if 
return 


end 


subroutine biterriiafill(data,nbits,insync,nerror) 


integer*4 data, nbits, nerror, rxdata, shindex, error 
integer*4 index, errcnt, bitcnt, i, itemp 
integer*4 nblck, btest, lshift, rshift, xor 
logical insync 
data errent, bitcnt / 0, 0 / : 
data rxdata, index / 0, 1 / 
data nblck /100/ 


rxdata = rshift(rxdata,nbits) 
call mvbits(data,0,nbits,rxdata,11-nbits) 


shindex = lshift(index,2) 

itemp = xor(index,shindex) 

index = rshift(index,nbits) 

call mvbits(itemp,2,nbits,index,11-nbits) 


nerror = 0 

error = xor(index,rxdata) 

if Cerror .ne. 0) then 7 
do 11 i =li-nbits, 10 


if (btest(error,i).ne.0) nerror = nerror +1 
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11 continue 
end if 


errent = errcnt + nerror 
bitcnt = bitcnt + nbits 


if (insync) then 
if (bitcnt .ge. nblck) then 
if (float(errcat) .gt. .35 * float(bitcnt)) then 


insync = .false. 
nbick=100 
else 
write(6,*) ’errors in ’,nblck, ’ bits =’, errcnt 
nblck=3000 
c to span fades 
errcent = 0 
bitcnt = 0 
end if 
end if 


end if 


if (.not. insync) then 
if (bitcnt .ge. 20) then 
if (float(errcnt) .le. .2 * float(bitcnt)) then 
insyne = .true. 
write(6,*) ~>insync = true’ 
else 
if (bitcnt.eq.20) write(6,*) ’not yet in sync’ 
index = rxdata 
if (index .eq. 0) index = 1 


end if 
errent = 0 
bitcnt = 0 
end if 
end if 
return 


end 
CR ROR ROR kk kok RR mig IO ROO OK ROR koko kok kk gaioiok ii ggoiok ak ak ak ae ak 
function btest(i1,i2) 


C *** this function returns 1 if bit position a2 of al is 1 and 
C #** returns 0 if bit position a2 of ai is 0 
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integer*4 t1,t2,btest,i1,i2 
tisil 

t2=1i2 

ti = rshift(t1,t2) 

ti = and(ti,1) 

btest = ti 


return 
end 


CoO IC OR OG i oR om igi ie gigi io i IGK kkuk kok kok ak i kok ak ak ak ae 


subroutine mvbits(isource,ibitfrom,length,idest,ibitto) 


integer*4 isource, ibitfrom, length, idest, ibitto, n 
integer*4 k, ks, btest, xor, lshift 


do 12 n= 1, length 


k = ibitto+n- 1 
ks = lshift(1,k) 


if (btest(isource, (ibitfromtn-1)).ne.0) then . 
Cc **sat** 
if (btest(idest,k).eq.0) 
1 idest = xor(idest,ks) - 
c **clear** 
else 
if (btest(idest,k).ne.0) 
1 idest = xor(idest,ks) 
end if 
12 continue 
return 
end 
a a ee ae 
c soft decision decoder for golay(24,12) 
uses chase ii algorithm with table look-up decoder . 


subroutine soft 





13 


anaana 


common/b1k5/idata 
common/b1k7/conf ,isym,kdata,ierr 
integer*4 idata(24) ,kdata(24) ,ierr(5) 
integer*4 isym(2047) 
integer*4 nb(4) ,mask(4,24> ,jb(16,24) ,ib(24) ,nerr 
real conf (24) 


nerr=i 


w is analog weight of most likely error pattern 
set w initially to a very large value 
w=1000000.0 


identify 4 bits with lowest confidence values 
call low4(nb) 
wgt4=0. 


gen 4 masks identifying 4 bits 
call mgen(mask ,nb) 


gen 2**4=16 patterns 
call patgen(mask, jb) 


test each of 16 patterns for the SOFT decoder 
do 19 i=1,16 


clear ierr 

do 13 j=1,5 
ierr(j)=0 

continue 


add data and test error pattern, store in kdata 
do 14 j=1,24 
call exor(kdata(j) ,idata(j),jb(i,j)) 
continue 


binary decoder 

accepts input data in kdata 

returns data uncorrected 

returns no. of errors in ierr(5) 

returns location of errors in ierr(1) to ierr(4) 
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call golay 


c test if decoded word was error free 
if (ierr(5).eq.0) then 
test word was error free 


c accept test word as best estimate 
w=0.0 
do 15 j=1,24 

idata(j)=kdata(j) 
15 continue 
return 
else 

test word contained errors 

c get final error pattern, calculate its analog wgt 
wgt=0.0 
n=ierr(5) 

c add detected errors to test error pattern 
do 16 j=1,n 


ix=ierr(j) 
call exor(jb(i,ix),jb(i,ix) ,nerr) 


16 continue 
do 17 j=1,24 
if(jb(i,j).eq.1)wgt=wgt+conf (j) 
17 continue 


compare wgt and w 
save error pattern with lowest analog weight 


if (wgt.1t.w)then 

c store new value 
save error pattern 
wewgt 
do 18 j=1,24 

ib(j)=jb(i,j) 
18 continue 
end if 
end if 
19 continue 


c decode data with most likely error pattern 


do 21 i=1,24 
call exor(idata(i),idata(i),ib(i)) 
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21 continue 





return 
end 
Ho pode eaboGtceese desl a pcceucesaseocgntSansatan wae ee asee enema = oe ae eeeeseeees 
C hard decision decoder for golay(24,12) 
c uses chase ii algorithm with table look-up decoder 
subroutine hard 
integer*4 idata(24) ,kdata(24) ,ierr(5) 
integer*4 isym(2047) 
integer*4 nb(4) ,mask(4,24) ,jb(16,24) ,ib(24) ,nerr 
real conf (24) 
common/b1k5/idata 
common/b1k7/conf ,isym,kdata,ierr 
nerr=1 
c wis analog weight of most likely error pattern 
c set w initially to a very large value 
w=1000000.0 
c identify 4 bits with lowest confidence values 
call low4(nb) 
wgt4=0. 
c gen 4 masks identifying 4 bits 
call mgen(mask ,nb) 
c gen 2**4=16 patterns 
call patgen(mask, jb) 
c test 1 pattern --> hard 
do 28 i=1,1 
c clear ierr 
do 22 j=1,5 
jerr(j)=0 
22 continue 
c add data and test error pattern, store in kdata 
do 23 j=1,24 
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call exor(kdata(j),idata(j) ,ib(i,j)) 


23 continue 


aaaonaana 


binary decoder 

accepts input data in kdata 

returns data uncorrected 

returns no. of errors in ierr(5) 

returns location of errors in ierr(1i) to ierr(4) 


call golay 


c test if decoded word was error free 


if 


(ierr(5).eq.0) then 


c test word was error free 
c accept test word as best estimate 
w=0.0 
do 24 j=1,24 
idata(j)=kdata(j) 
24 continue 
return 
else 
c test word contained errors 
c get final error pattern, calculate its analog wgt 
wgt=0.0 
n=ierr(5) . 
Cc add detected errors to test error pattern 
do 25 j=i,n 5 
ix=ierr(j) 
call exor(jb(i,ix) ,jb(i,ix) ,nerr) 
25 continue 
do 26 j=1,24 
if (jb(i,j) .eq. 1) wgt=wet+conf (j) 
26 continue 
c compare wgt and w 
c save error pattern with lowest analog weight 





if (wgt.1t.w)then 
store new value 
save error pattern . 
wewgt 
do 27 j=1,24 











ib(j)=jb(i,j) 
27 continue 
end if 
end if 
28 continue 


c decode data with most likely error pattern 
do 29 i1=1,24 
call exor(idata(i),idata(i),ib(i)) 


29 continue 
return 
end 

c 


subroutine low4(nb) 


C to identify in nb the location of 4 lowest conf values 
dimension d(4) 
integer*4 kdata(24),ierr(5) ,nb(4) 
integer*4 isym(2047) 
real conf (24) 
common/b1k7/conf ,isym,kdata,iei > 


c assume first 4 are lowest 


do 31 i=1,4 
d(i)=conf (i) 
nb(i)=i 
31 continue 


c compare remaining values 
do 33 125,24 
a=conf (i) 
n=i 
j=l 
if (j.gt.4) go to 33 
if (a.1t.d(j)) then 
ais < d(j) 
replace d(j) with a 
save d(j) to compare with other values in d 
save nb(j) 
b=d/j) 
nx=nb(j) 


a an a 








33 


34 


35 





replace 


d(j)=a 
nb(j)=n 
rename ihe saved vaiues 
a=b 
n=nx 
j=1 
end if 
incr j 
jejti 
continue 
return 
end 


mask generator for soft decision golay decoder 
generate 4 masks with 1 bit in location nb(i) 


subroutine mgen(mask,nb) 


dimension mask(4,24) ,nb(4) 
integer*4 mask ,nb 


do 35 i=1,4 
do 34 j=1,24 
mask(i,j)=0 
continue 
n=nb(i) 
mask(i,n)=1 
continue 


return 
end 


_—“—Ke wn Oe OO HOw Oe we wwe mw ew mw Me ee we eww mw em me eww mem we we ew eee ee ee ewe 


generate 16 patterns based on 4 bits 
subroutine patgen(mask, jb) 
integer*4 mask(4,24) ,jb(16,24) 


do 36 1=1,24 

















36 


37 
56 





jb(1,i)=0 
jb(2,i)=mask(1,i) 
jb(3,i)=mask(2,i) 
jb(4,i)=mask(3,i) 
jb(5,i)=mask(4,i) 
jo(6,i)=jb(2,i)+mask(2,i) 
jo(7,i)=jb(2,1)+mask(3,i) 
jo(8,i)=jb(2,i)+mask(4,i) 
jb(9,i)=jb(3,i)+mask(3,i) 
jb(10,i)=jb(3,i)+mask(4,i) 
jo(11,i)=jb(4,i)+mask(4,i) 
jb(12,i)=jb(6,i)+mask(3,i) 
jb(13,i)=jb(6,i1)+mask(4,i) 
jb(14,i)=jb(7,i)+mask(4,i) 
jb(15,i)=j5(9,i)+mask(4, i) 
jb(16,1)=jb(12,i)+mask(4, i) 
continue 


do 56 j=1,16 
do 37 i=1,24 
jo(j,i)=and(jb(j,i),1) 
continue 
continue 


return 
end 


Ste ee em ee ww wr em wee mee mw em we mm em em em mm me mmm meme we em mm mm me eee ew mm em ew we em ee ee eee ee 


exclusive or of two integer*4s 
one bit per integer*4 


subroutine exor(i,j,k) 
integer*4 i,j,k 
isand(j+k,1) 


return 
end 


we mm em Me me em ee Hw Pee mw tet ew ew ew we www ow ew ew ew ew em em ew ew ew ee ew ew ew ew ee eee ew eee ee eee ee ee 


binary table look-up decoder for golay (24,12) code 











subroutine golay 


integer*4 idata(24),ierr(5) ,kdata(24) ,is(11) 
integer*4 isym(2C47) 

real conf (24) 

common/b1k5/idata 

common/b1k7/conf ,isym,kdata,ierr 


n=23 
ierr(5)=0 
Cc generate syndrome 
c write(6,*) ’generating golay syndrome’ 


call encode(kdata,is,n) 


c use syndrome to generate index to table 
ir=0 
do 38 j=i,11 
ir=ir*2+is(j) 


38 continue 


if(ir.ne.0)then 
look up error pattern 
possibility of up to 3 errors 


Cc recorded as three 5-bit symbols in a 16-bit word 
ix=isym(ir) 
C identify error locations and count errors 


do 39 j=1,3 
ierr(j)=and(ix,31) 
if(ierr(j).gt.0)ierr(5)=ierr(5)+1 


ix=1ix/32 
c if(ierr(j).gt.0)write(6,*)’ierr(j)=’ ,ierr(j) 
39 continue 
end if 
¢C test overall parity 
ir=ierr(5) 
do 41 j=1,24 
ir=ir+kdata(j) 
41 continue 
ir=and(ir,1) 
c test for even parity 


if (ir.ne.0) then 


parity check failed 
force overall parity 
increment error count 
ir=ierr(5)+1 

jierr (ir) =24 
jerr(5)=ir 


return 


golay(23,12) encoder, syndrome generator 


subroutine encode(idata,ip,n) 
integer*4 idata(12),ip(11),ix 


F(X) EKHRL LEX EROS KEK +X HRETXEHS + X41 
clear ip 
do 42 i=1,11 
ip(i)=0 
continue 


read inn bits 

do 43 i=1i,n 
call exor(ix,idata(i),ip(1)) 
ip(1)=ip(2) 
call exor(ip(2),ix,ip(3)) 
ip(3)=ip(4) 
call exor(ip(4),ix,ip(5)) 
call exor(ip(5),ix,ip(6)) 
call exor(ip(6),ix,ip(7)) 
ip(7)=ip(8) 
ip(8)=ip(9) 
ip(9)=ip(10) 
call exor(ip(10),ix,ip(11)) 
ip(11)=ix 

continue 


return 








subroutine bitgenii(nbits,ival,data) 


this routine generates a pseudorandom sequence of bits using an 11 bit 
shift register, exoring the Sth and 11th bits to generate the new ist bit 
before shifting. the length of the pseudorandom sequence is (2**11)-1 

or 2047 bits. 


aanan 


integer*4 data, shindex, nbits, ival, index, temp 
integer*4 lshift, rshift, xor 


if (ival .gt. 0 .and. ival .1t. 2048) index = ival 


data = 0 
c ensure that all bits in data are 0 


call mvbits(index,0,nbits,data,0) 
c move bits from index to data 


shindex = lshift (index, 2) 
c shindex is index shifted left 2 


temp = xor(index,shindex) 
c temp contains new bits for index 


index = rshift(index,nbits) 
shift index right by nbits 


fe) 


call mvbits(temp,2,nbits,index,11-nbits) 
c move new bits, temp to index 


return 
end 


c golay encoder 
subroutine golenc 


integer*4 idata(24) ,id(12),is(11) 
common/b1k5/idata 


n=12 





ce load 12 information bits 
c sum inf bits for overall parity 


do 44 j=1,12 
id(j)=idata(j) 
44 continue 


c encode using n-k type shift register 
call encod(id,is,n) 


c store parity bits in same array as inf bits 


do 45 j=1,11 
idata(j+12)=is(j) 
45 continue 


© generate over all parity bit(even parity) 
c store as 24th bit in array 
m=0 


do 46 j=1,23 
m=mt+idata(j) 
46 continue 


idata(24)=and(m, 1) 


return 
end 


¢ golay(23,12) encoder, syndrome generator 
C x114+x9+x7+x6+x5+x+1 


subroutine encod(idata,ip,n) 
integer*4 idata(12),ip(11),ix 


c clear ip 


do 47 i=1,11 
ip(i)=0 
47 continue 


c read inn bits 
do 48 i=1,n 
call exor(ix,idata(i),ip(1)) 
ip(1)=ip(2) 
call exor(ip(2),ix,ip(3)) 








ip(3)=ip(4) 

call exor(ip(4),ix,ip(5)) 
call exor(ip(5),ix,ip(6)) 
call exor(ip(6),ix,ip(7)) 
ip(7)=ip(8) 

ip(8)=ip(9) 

ip(9)=ip(10) 

call exor(ip(10),ix,ip(11)) 


ip(11)=ix 
48 continue 
return 
end 
Caseeesscb ies socccs ou ons cece se scki ng co eksccsetcsesee cess seeuSestsecsecsssusalle 
subroutine bitgeniifill(nbits,ival,data) 
c this routine generates a pseudorandom sequence of bits using an 11 bit 
c shift register, exoring the 9th and 1ith bits to generate the new 1st bit 
c before shifting. the length of the pseudorandom sequence is (2**11)-1 
c or 2047 bits. 


integer*4 data, shindex, nbits, ival, index, temp 
integer*4 lshift, rshift, xor 


if (ival .gt. 0 .and. ival .1t. 2048) index = ival 


data = 0 
c ensure that all bits in data are 0 


call mvbits(index,0,nbits,data,0) 
c move bits from index to jiata 


shindex = lshift (index, 2) 
c shindex is index shifted left 2 


temp = xor(index,shindex) 
c temp contains new bits for index 


index = rshift(index,nbits) 


c shift index right by nbits 
call mvbits(temp,2,nbits,index,1i-nbits) 
c move new bits, temp to index 








return 
end 


subroutine tx(dibit,txsamp) 


dimension txbuff(32), txsamp(3) 
integer*4 dibit, index, i 

rea’ txbuff, txsamp 

data txbuff, index / 32*0., 1 / 


call modulate(dibit,index,txbuff) 
c modulate one baud or symbol 


c add symbol to transmit buffer 
c note : each symbol is filtered by the modulator to be wider than 
¢ one baud so as to restrict the bandwidth to about 4000 Hz for 
c reduction of ACI (adjacent channel interference), so that the resulting 
c baud pulses overlap and must be added together 
do 49 i= 1, 3 
txsamp(i) = txbuff (index) 
txbuff (index) = 0 
index = mod(index,32) + 1 
49 continue 
return 
end 
CO ee ee meme mew we ew we mem em em ew ee ee eee ee ee ee 


subroutine modulate(mbit, index,buffer) 


integer*4 mbit, phasec, cphase, tphase 
integer*4 point, wlength, i, index, modangle 
real s, pi, buffer, w 

dimension phasec(4), w(29), buffer(*) 


data phasec / 225,135,315,45 / 
data wlength, modangle, cphase, pi / 29, 0, 0, 3.141593 / 


c this is the transmit window to shape the output spectrum to about 
c 4000 Hz wide 











data w / -.005696, .062601, .028191,-.021020,-.060997, 
- .035666, .049659, . 106424, .043728, - . 114226 ,- .203894, 
-.049151, .358925, .805832, .999756, .805832, .358925, 

- .049151,-.203894,-.114226, .043728, .106424, .049659, 
- .035666 , - .060997 ,- .021020, .028191, .062601 ,- .005696/ 


* % & & 


modangle = phasec(mbit + 1) 
modulation phase change angle 


cphase = mod(cphase + 270, 360) 
continuous unmodulated carrier phase 


tphase = mod(cphase + modangle, 360) 
actual transmit phase 


point = index 
take the sine wave of the actual transmit phase and window it 
do 51 i = 1, wlength 
s = .364 * w(i) * sin(pi / 180. * (tphase + 90 * i)) 
buffer(point) = buffer(point) + 8 
point = mod(point, 32) + 1 


continue 


return 
end 


subroutine rcv(samples,dbits,xdisp,ydisp,confh,confl) 


, purpose: this routine demodulates an 8000 b/s 4 phase dpsk modem signal. 
it is called once to demodulate each baud. 


inputs: sigin - contains samples of the signal to be demodulated. 
sigin(1) is first in time. 


outputs: dbits - an integer which contains the 2 demodulated bits. 
The data is stored in the least significant bits; 
the lsb is first in time. 


brief description of subroutine functions: 


equalize - calculates the output of the adaptive equalizer. 
decode - uses the equalizer output to determine the data bits. 


80 


a phase reference is calculated for internal use. 


integer*4 dbits, i 
real confh,confl 

reai signal, x, y 

real xdisp,ydisp 

real signal(48), samples(3) 
real xcoef (48) ,ycoef (48) 


common /adequ/ xcoef, ycoef 





c compromise equalizer coefficients for remeving the intersymbol 
c interference introduced by the transmitter windowing : 
data xcoef / 0.00047277, 0.00018969, 0.00152612,-0.00026510, 


1 0.00288611,-0. 
-0.00299245 ,-0. 
.00073868, 0. 
.01571514,-0. 
.03999998 ,-0. 
.14213182, 0. 
.06119568, 0. 
.01294118, 0. 
.01742128, 0. 
-00697080,-0. 
.00324149,-0. 


00077746, 0. 
00031443,-0. 
00575836, 0. 
00732426, 0. 
00832545 ,-0. 
01429809, 0. 
01810114,-0. 
00074462, 0. 
00060S89,-0. 
00119318, 0. 
00039312, 0. 


00057433 ,-0. 
01333420, 0. 
00582585, 0. 
00609350,-0. 
10446436, 0. 
38911167, 0. 
00429513, 0. 
00006572, 0. 
00858223, 0. 
00344639 ,-0. 
00125349,-0. 


00176436, 
00245006, 
00128170, 
01291209, 
00526107, 
03028025, 
00713205, 
00038059, 
00059077, 
00073849, 
00045128/ 





data ycoef / -0.00016221, 0.00070496 ,-0.00036193,-0.00072445, 
-0.00059857 ,-0.00178603 ,-0.00014448 ,-0.00418214, 
0.00070030, 0.00100549, 0.00305764, 0.00433583, 
0.00079199, 0.00792035 ,-0.00374299, 0.01238339, 
-0.00632491 ,-0.01299071, 0.00042162,-0.03951978, 
0.01424530 ,-0.02844719, 0.02507262, 0.05758527, 
0.02317356, 0.35042825, 0.01940719,-0.24283291, 
~0.00569436 ,-0.01806669 , -0.01134070 ,-0.00663402, 
~0.00711616,-0.02175912,-0.0007£.578,-0.01645218, 
0.00256463 ,-0.00212019, 0.00016785, 0.02115092, 
-0.00076241, 0.00839603 ,-0.00143184,-0.00148724, 
-0.00056496 , -0 .00391606 , -0.00031074,-0.00173826/ 





fe eo ee ee ee ee ee ee ee ee 


data signal / 48*0.0 / 


¢ signal buffer is as long as the equalizer that it will multiply against 
c or 48 samples long 


do 52 i= 1, 48 
shift signal buffer to accomodate new samples 
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signal(i) = signal(i+3) 
52 continue 


do 53 i = 46, 48 


c put new samples into signal buffer 
signal(i) = samples(i-45) 
53 continue 


call equalize(signal,x,y) 
call decode(x,y,dbits,xdisp,ydisp,confh,conf1l) 


return 
end 


subroutine equalize(signal,x,y) 


¢ purpose: this routine generates the output of the equalizer; i.e., it 
c ‘culates the output value of each fir equalizer filter. 


c input: signal - a real array which contains the received signal samples. 
c signal(i1) is first in time. 


c outputs: x - the x-coordinate of the equalizer output. 
y - the y-coordinate of the equalizer output. 


integer*4 i 
real x, y 

real xcoef (48) ,ycoef (48) 
real signal (48) 


common /adequ/ xcoef, ycoef 


y =0. 
x= 0. 


do 54 i = 1,48 
x = x + signal(i) * xcoef(i) 
y = y + signal(i) * ycoef(i) 
54 continue 


X=2.*X 
y=2.*y 


io 8) 
we 





return 
end 

Qoetese cece eeeece ct eotee ce ceese tee lee Ss eae nS eb Seco ct ease se ceecees 
subroutine decode(x,y,dibit,xdisp,ydisp,confh,conf1) 

s purpose: this routine determines the data bits given the x,y coordinates 
from the equalizer. 
c inputs: x - the x-coordinate of the demodulated baud. 
c y - the y-coordinate of the demodulated baud. 





c outputs: dibit - the 2 data bits for the current baud; lsb is first in time. 


integer*4 dibit 
real mag, phaser, x, y, qang, preang, phasech, arctan 
real angle, decang, angdisp, xdisp, ydisp, error 
real angtab 
dimension angtab(4) 
real predecang,diffang 
real confh,confl 
real desmag 


data angtab / -.75, -.25, .25, .75 / 
data preang / 0. / 


angle = arctan(y,x) 
¢ convert from rectangular to polar coordinates 


decang = angle - phaser 
c decang is the angle to be decoded 


call adjust (decang) 
c calculate dibit: 


c preang=qang 
c previously quantized angle 


c qang = angtab(int(2.*decang+2.)+1) 
c quantized angle 


c phasech = qang - preang 





¢ quantized phase change 


Cc if (phasech.1t.0.0)phasech=phasech+2. 
¢ if (phasech.eq.0.0)dibit=2 

c phase change of 0 

c if (phasech.eq.0.5)dibit=0 

c phase change of 90 

c if (phasech.eq.1.0)dibit=1 

c phase change of 180,or -180 

c if (phasech.eq.1.5)dibit=3 

c phase change of -90 


diffang=decang-predecang 
c unquantized phase change 


call adjust (diffang) 


predecang=decang 
c previous decoded angle update 


c compute confidence values : 
if (diffang.1t.-.75)then 
dibit=1 
confl=(diffang+1.25)*2. 
confh=-(diffang+.75)*2. 
else if(diffang.1t.-.25)then 
dibit=3 
confl=-(diffang+.25)*2. 
confh=(diffang+.75)*2. 
else if(diffang.1t.0.25)then 
dibit=2 
confl=(diffang+.25)*2. 
confh=~(diffang-.25)*2 
else if(diffang.1t.0.75)then 
dibit=0 
confl=-(diffang-.75)*2. 
confh=(diffang-.25)*2. 
else 
dibit=1 
confl=(diffang-.75)*2. 
confh=-(diffang-1.25)*2. 
end if 


c calculate x and y values for constellation display: 
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mag = sqrt(x*x + y*y) 
angdisp = (diffang+.25) 

call adjust (angdisp) 

angdisp = angdisp * 3.141593 


xdisp = mag * cos(angdisp) 
ydisp = mag * sin(angdisp) 


desmag=.707/2 
c scale the confidence by the radius (lower confidence in fade) 
if (mag .1t .desmag) then 
confl=confl+*mag/desmag 
confh=confh*mag/desmag 
end if 


c calculate error angle : 


c error = decang - gang 
c call adjust (error) 


c update phase reference for next baud: 


phaser = phaser - .5 


c the .5 comes from .25 (times 180 degrees) 
c of carrier advance per baud plus .25 
c (cr 45 degrees) thus shifting the phase computations 
c from 135,45,-45,-135 to the 90,0,-90,-180 domain 
call adjust (phaser) 
return 
end 
Cc eS ee Ee eS ee ee a eee ee ee ee ee ee ee eee eee ee ee eee ee eee ee ewe a ee ee 


function arctan(y,x) 


c 4 quadrant arctangent -1 <= arctan < 1 
c -1 corresponds to -180 degrees 


real arctan, y, x 


arctan = 0.0 
if (x .ne. 0.0 .or. y .ne. 0.0) arctan = atan2(y,x)/3.141593 


return 
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subroutine adjust (value) 

¢ routine to ensure value is kept in range -1 to 1. 
real value 
value = amod(value,2.) 


if (value .1t. -1.) value = value + 2. 
if (value .ge. 1.) value = value - 2. 


return 
end 


c noise=signal/noise in dB 
c mode = 0 for non-fading, 1 for fading 


real samples(3) 
real gaus(256) 
c gaussian distribution of noise table 
real noise,noisef 
real sv(120) ,nv(120) 
real sdb,ndb 
integer*4 sinn,nin,n,sf 
integer*4 fade,mode,fadecnt 
integer*4 iseed 
integer*4 ms10,ms20,ms40,ms100,ms200 


if (ifirst.eq.0) then 
rgn = (10.**(-noise/20.)) 
c compute actual noise scaler from dB 


rgnf = (10.+**(-noisef/20.)) 


c compute noisef scaler from dB 
ifirst=1 
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test=rand(iseed) 


fade = 0 

fadecnt = 0 

msiO = 0 

ms20 = 0 

ms40 = 0 

msi00 = 0 

ms200 = 0 
end if 
a a ee ae a a a ee a 
do 55 m=1,3 


c measure the RMS of past 120 samples of signal 
s = samples(m) 
sinn = mod(sinn+t1,120) 
sdb = (sdb**2.)-sv(mod(sinn+1,120)+1) 
sv(sinn+1) = (s**2.)/120. 
sdb = sdb + sv(sinn+1) 
if (sdb.1t.0.) sdb = 0. 
sdb = sqrt(sdb) 


c measure the RMS of past 120 samples of noise (direct from table) 
nin = mod(nin+1,120) 
ndb = (ndb**2.)-nv(mod(nin+1,120)+1) 
k = (255.*rand(0)) + 1 
nv(nint+i) = ((gaus(k)/25295.)**2.)/120. 
ndb = ndb + nv(nin+1) 
if (ndb.1t.0.) ndb = 0. 
ndb = sqrt(ndb) 


if (mode.eq.1)then 
if (fadecnt.eq.0)then 
if (rand(0) .gt..94056) then 
fade=1 
test=rand(0) 
if (test.1t.0.8)then 
fadecnt=120 
msi0 = msi0 + 1 
else if(test.1t.0.9)then 
fadecnt=240 
ms20 = ms20 + 1 
else if(test.1t.0.95)then 
fadecnt=480 
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ms40 = ms40 + 1 
else if(test.1t.0.99)then 
fadecnt=1200 
msi00 = msi00 + 1 


else 
fadecnt=2400 , 
ms200 = ms200 + 1 


end if 
else 
fade=0 
fadecnt=120 
end if 
else 
fadecnt=fadecnt-1 
end if 


if ((n.eq.sf) .and. (mk.eq. 1680) .and. (m.eq.3))then 
write(6,*)’ number of 10ms fades in ’,sf, 
> ,ms10 
write(6,*)’ number of 20ms fades in ’,sf, 
? ,ms20 
write(6,*)’ number of 40ms fades in ’,sf, 
> ,ms40 
write(6,*)’ number of 100ms fades in ’,sf, 
*’ runs = ’,msi00 
write(6,*)’ number of 200ms fades in ’,sf, 
* ,ms200 
end if 
end if 


*? runs 


*’ runs 


*? runs 


*? runs 


if ((mode.eq.0).or. ((mode.eq.1).and. (fade.eq.0))) then 
c not in fade 
if (rand(0).ge.0.5) then 
if (ndb.ne.0.) samples(m) = samples(m) + 


1 (gaus(k)/25295.)*(sdb/ndb) *(rgn) 
else 
if (ndb.ne.0.) samples(m) = samples(m) - 
1 (gaus (k) /25295.)*(sdb/ndb) *(rgn) 
end if 
end if 


if ((mode.eq.1).and. (fade.eq.1))then 
c we are in fade 
if (rand(0).ge.0.5) then 
c make positive noise 








if (ndb.ne.0.) samples(m) = samples(m) + 
1 (gaus (k) /25295.)*(sdb/ndb)*(rgnf) 
else 
c make negative noise 
if (ndb.ne.O.) samples(m) = samples(m) - 
1 (gaus(k)/25295.)*(sdb/ndb) *(rgnf) 
end if 
end if 


c scale signal down in fade : 
if ((mode.eq.1).and. (fade.eq.1))samples(m)=samples(m)/10. 


55 continue 


return 
end 


COR RK Rk kkk OR kok igi RO Rk kok gk Rk GR Rk kkk kkk kok kok ak ak kkk kok kak kok kk koko kak kkk ok ak 


subroutine chsim(samples ,noise,noisef ,mode,iseed,gaus,bprob,bwidth) 


c noise=signal/noise in dB 
c mode = O for non-fading, 1 for fading 


real samples(3) 
real gaus(256) 
c**ee* gaussian distribution of noise table 
real noise,noisef ,bprob 
real sv(120) ,nv(120) 
real sdb,ndb,rtest 
integer*4 sinn,nin 
integer*4 fade,mode,fadecnt 
integer*4 iseed,bwidth 


if (ifirst.eq.0) then 
rgn = (10.**(-noise/20.)) 
c**e*e* Compute actual noise scaler from dB 


rgnf = (10.**(-noisef/20.)) 
c###** compute noisef scaler from dB 
ifirst=1 
test=rand(iseed) 
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fade = 0 


fadecnt = 0 

rtest = 1. - (1.-bprob)/120 
end if 
prepa een ankoen cde obese oases sess nte oct ata se ao she etc 
do 55 m=1,3 


c***#* measure the RMS of past 120 samples of signal 
s = samples(m) 
sinn = mod(sinn+1,120) 
sdb = (sdb**2.)-sv(mod(sinn+1,120)+1) 
sv(sinnt1) = (s**2.)/120. 
sdb = sdb + sv(sinn+i) 
if (sdb.1t.0.) sdb = 0. 
sdb = sqrt(sdb) 


c**#** measure the RMS of past 120 samples of noise (direct from table) 
nin = mod(nin+1,120) 
ndb = (ndb**2.)-nv(mod(nin+1,120)+1) 
k = (255.*rand(0)) + 1 
nv(nin+1) = ((gaus(k)/25295.)**2.)/120. 
ndb = ndb + nv(nin+1) 
if (ndb.1t.0.) ndb = 0. 
ndb = sqrt(ndb) 


if (mode.eq.1)then 
if (fadecnt.eq.0)then 
if (rand(0) .gt.rtest)then 
fade=1 
fadecnt = bwidth 
else 
fade=0 
fadecnt=0 
end if 
else 
fadecnt=fadecnt-1 
end if 
end if 


c*x**ke* Not in fade 


if ((mode.eq.0) .or. ((mode.eq.1).and. (fade.eq.0))) then 
if (rand(0).ge.0.5) then 
if (ndb.ne.0.) samples(m) = samples(m) + 
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1 (gaus(k)/25295.)*(sdb/ndb) *(rgn) 


else 
if (ndb.ne.0.) samples(m) = samples(m) - 
i (gaus(k) /25295.)*(sdb/ndb) *(rgn) 
end if 
end if 


c***k* we are in fade 
if ((mode.eq.i) . and. (fade.eq.1))then 
c*x**e* make positive noise 
if (rand(0).ge.0.5) then 
if (ndb.ne.0.) samples(m) = samples(m) + 
1 (gaus(k) /25295.)*(sdb/ndb) *(rgnf) 
cxxke* make negative noise 
else 
if (ndb.ne.0.) samples(m) = samples(m) - 
1 (gaus(k)/25295.)*(sdb/ndb) *(rgnf) 
end if 
end if 


c***** scale signal down in fade : 
if ((mode.eq.1).and. (fade.eq.1))samples(m)=samples(m)/200. 


55 continue 


return 
end 


COR ROR doko kom mk kokok kkk ok 


c encodeham 

Cc 

c FUNCTION 

c This subroutine calculates the parity bits necessary 
c to form the codeword. 

c 

c 

c SYNOPSIS 

c encodeham(codelengthi,codelength2,hmatrix, 
c paritybit ,codeword) 

c 

¢ formal 

c 

c data 1/0 

c name type type function 
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codelengthi int i number of data bits (63) 
codelength2 int i number of information bits (57) 
hmatrix int i vector to encode an decode by 
paritybit int o overall parity bit 
codeword int o encoded stream (paritybits at end) 
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DESCRIPTION 


This subroutine is part of a set of subroutines which perform 

a Generalized Hamming Code. As you know, Hamming codes are perfect 
codes and can only detect and correct one error. We added an overall 
parity checkbit, which allows us to detect 2 errors. When 2 errors 
are detected, (in subroutine decodeham.f) no correction attempt is 
made. This would most likely result in more errors. Instead, a flag 
is sent to the calling program notifying it of multiple errors so 

that smoothing may be attempted. The Hamming codes presently supported 
by the routines are (63,57), (31,26), (15,11), and shortened variations 
thereof. It could be made even more general by making minor modifications 
to the dectobin.f subroutine. This routine at present will calculate 

a maximum of 6 bits. 


Hamming routines consist of the following files: 


Matrixgen - generates the hmatrix and sydrometable. 

dectobin - does a simple decimal to binary conversion. 
encodeham - generates the codeword and overall paritybit. 
decodeham - recovers infobits, checks for errors, corrects i 
error, and sends out flag for smoothing. 


This subroutine performs the Hamming encode function. 
It will calculate the necessary parity bits, depending on which code 
is requested, and will add the overall parity bit to the end of the 
codeword generated. 


aoaoanqaoanoanqananoanaoaonanaoaqaoqenandneaanqnaanaaoaoaqaaoaqoaoa0o0aasa 


CSS SSS SSS SSSSSS SSP SSSSSP SSS SSSSSSSSSSSSssssseqssssseresseessseesrssssserrq==e2=== 

c 

c REFERENCES 

v ° 
c¢ Lin and Costello : Error Control Coding 

c Berlekamp : Algebraic Coding Theory 
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Cc 
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c 

subroutine encodeham(codelengthi ,codelength2,hmatrix, 
ae paritybit ,codeword) 

c 


integer*#4 codelength1 ,codelength2,paritybit 

integer*4 codeword(codelength1) ,hmatrix(codelength1) 
integer*4 parityflag,temp1,temp2,i,temp3 

c 

paritybit=0 

parityflag=0 

c *** parityflag = 0 if not using the extra parity bit. 
temp1=codelengthi-codelength2 


temp2=0 

temp3=0 

c 

c First generate the parity bits for the Hamming codeword. This is 
c relatively straightforward. hmatrix was generated in matrixgen.f, 
c which is called as part of the Hamming initialization routines. 

c 


do 10 i=1, codelength2 
if (codeword(i).ne.0) temp2= xor(temp2,hmatrix(i)) 
10 continue 


since the hmatrix is stored in a packed decimal format, the parity 
bits must be unpacked and appended to the end of the bitsteam. 
after this routine you will have the complete codeword. 


aaaaa 


call dectobin(temp1,temp2, codeword (codelength2+1)) 


c 

c Now I check to see if the parityflag is set, indicating the user 
c requests an overall parity bit be generated. Normally this will 
c be the case. 

¢c 

temp2=0 


if (parityflag.eq.1)then 
do 20 i=1,codelengthi 
temp2 = xor(temp2,codeword(i)) 
if (codeword(i).ne.0) temp3=temp3+1 


20 continue 
paritybit=temp2 

end if 

c 


return 





ROUTINE 
decodeham 


FUNCTION : 
This subroutine decodes the bitstream generated by 

encodeham. It will correct a single error, and detect 2 

errors. 


SYNOPSIS 
subroutine decodeham(codelength1 ,codelength2,hmatrix, 
syndrometable,paritybit ,codeword,myerror) 


formal 


data 1/0 


we ww em we em mw mew em mem wee ee mew we em ee ewe ww ww wwe we mw wwe Hew ww ww ee em ee wee eee 


codelengthi int i number of data bits 
codelength2 int i number of information bits 
hmatrix int i vector to encode an decode by 
syndrometable int i errormasks used to correct single 
errors 
paritybit int i overall parity bit 
codeword int i/o encoded/decoded stream 
myerror log o flag for 2 error detect 
synflag int o value 0 or 1, 1 if syndrome .ne. 0 


Aanonaeadqaaaoeanaqonianoeananonanoanoenoeanqneaeaqaaqaonianqneaeaqaeniaoaaa 


Pe ee ee ee ee ee ee ee ee eee 
SPSS SSSSSS SSS SSS SS SSS SSSS SSS SSS SSS SSS SSSSS SSS SSeS SsrSsSsSS Sse essere sezsssescsrreza= 


a 


DESCRIPTION 


This subroutine is part of a set of subroutines which perform 

a Generalized Hamming Code. As you know, Hamming codes are perfect 
codes and can only detect and correct one error. We added an overall 
parity checkbit, which allows us to detect 2 errors. When 2 errors 
are detected, (in subroutine decodeham.f) no correction attempt is 

made. This would most likely result in more errors. Instead, a flag 
is sent to the calling program notifying it of multiple errors so 

that smoothing may be attempted. The Hamming codes presently supported 


aaanonaoaqaaanada 
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by the routines are (63,57), (31,26), (15,11), and shortened variations 
thereof. It could be made even more general by making minor modifications 
to the dectobin.f subroutine. This routine at present will calculate 

a maximum of 6 bits. 


Hamming routines consist of the following files: 


Matrixgen - generates the hmatrix and sydrometable. 

dectobin - does a simple decimal to binary conversion. 

encodeham - generates the codeword and overall paritybit. 

decodeham - recovers infobits, checks for errors, corrects 1 , 
error, and sends out flag for smoothing. 


This subroutine, decodeham, is resporsible for checking for errors, 
correcting the error if there is only one, and sending a smoothing fleg 
to the calling routine if there is more than one. 


oadndnoaqaeqnqnenqonaqeqoanqaoaeaqaaqoqonoaoeaana 


Set et i ee ww ow a a a we ee 
SS SS SSS SSS 2 SS SSS SSS SSS SS SS SS SS SSS SSS SSS SSSR SS SSSSSSSHS 22 5o EES 2qEE EES eESEsz5rrrEz==zm= 


a 


REFERENCES 


Lin and Costello : Error Control Coding 
Berlekamp : Algebraic Coding Theory 


angqaaanaa 


COI OO OO OR OI GK OR OR RO OO iI i RIK aK kok a 
c 
subroutine decodeham(codelengthi ,codelength2,paritybit ,myerror, 
1 hmatrix,syndrometable,codeword) 


integer*4 codelength1 ,codelength2,hmatrix(codelength1) 

integer*4 syndrometable(codelength1) ,paritybit 
integer*4 codeword(codelength1) ,myerror 

integer*4 errorflag,parityflag 

integer*4 synflag, i, j, temp3 


myerror=0 
errorflag=0 
parityflag=0 
*#e parity flag = 0 if not using the extra parity bit 


and compares it with the overall paritybit sent. If they are not the 
same that means there is at least one error. If, later on in the routine, 
the syndrome check indicates that there is an error and the parity is 





c 
c 
c This part of the routine checks the overall parity of the codeword 
c 
c 
c 


correct in this part of the routine, that indicates ther> are two errors. 
One of the weaknesses of this method is that there is no way of knowing 
if we have 3,5,7,... errors. We always smooth if there are 2,4,6,... 
errors. 


aaaana 


if (parityflag.eq.1) then 
synflag=0 
do 10 i=1,codelengthi 
synflag= xor(synflag,codeword(i)) 


10 continue 
if (paritybit.ne.synflag)errorflag=errcerflag+i 
end if 
c 
c This part of the routine generates the syndrome. The syndrome will 
c equal zero if there are no errors. synflag accumulates the syndrome 
c and is used as the offset in the syndrome table, which tells the 
c routine which bit is in error. 
Cc 


synflag=0 
temp3=0 
do 30 i=1,codelengthi 
if (codeword(i).ne.0)synflag = xor(synflag,hmatrix(i)) 
if (codeword(i).ne.0)temp3 = temp3 + 1 
30 continue 


c 
c **k* Check to see if the parityflag is set and if it is then check 
c to see if the parity bit was in error. 

c If the parityflag was set and there was an error in the syndrome, 
c the errorflag should equal 1. 

c If it doesn’t, then there are more errors than can be corrected 

c and the infobits are passed on unchanged. 


c 
if (synflag.ne.0)then 
if((errorflag.ne.1) .and. (parityflag.eq.1))then 
myerror = i 
go to 20 
end if 
j=syndrometable(synflag) 
codeword (j)=xor(codeword(j),1) 
end if 
c 
*x** If the syndrome is equal to zero and the errorflag is set 
(not likely, but must be checked) then more than one error has 
occured, but it cannot be corrected, so I pass on the infobits 
the same as if there were no errors. 


qa a0 0 
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20 continue 


r 


eturn 


end 


Onanqanaanananaanaanianianaaaaqaqaaoana 


aanoaqaqaaognoananoananaqenoaqaagqoea0oandna 





ROUTINE 
dectobin 


FUNCTION 
This subroutine converts decimal numbers into a 
binary output vector. 


SYNOPSIS 
dectobin(vectorsize, decinteger, binaryvector) 


formal 


data 1/0 
name type type function 
vectorsize int i output vector length 
decinteger int i decimal number (< 2732-1) 
binaryvector int o vector containing binary number 


DESCRIPTION 


This subroutine is part of a set of subroutines which perform 

a Generalized Hamming Code. As you know, Hamming codes are perfect 
codes and can only detect and correct one error. We added an overall 
parity checkbit, which allows us to detect 2 errors. When 2 errors 
are detected, (in subroutine decodeham.f) no correction attempt is 

made. This would most likely result in more errors. Instead, a flag 
is sent to the calling program notifying it of multiple errors so 

that smoothing may be attempted. The Hamming codes presently supported 

by the routines are (63,57), (31,2°), (15,11), and shortened variations 

thereof. It could be made even more general by making minor modifications 

to the dectobin.f subroutine. This routine at present will calculate 

a maximum of 6 bits. 


Hamming routines consist of the following files: 


matrixgen - generates the hmatrix and sydrometable. 


oT 





dectobin - does a simpie decimal to binary conversion. 
encodeham - generates the codeword and overall paritybit. 
decodeham - recovers infobits, checks for errcers, corrects 1 
error, and sends out flag for smoothing. 


This routine is used by encodeham to convert the packed decinteger 
into the Hamming paritybits. 


aanaaaanoaanda 


CAAA OI IGOR OIG ICSC OO IO IO IO GOIO IG ROI ISI II IIE 
c 

c REFERENCES 

c 

c Lin and Costello : Error Control Coding 

c Berlekamp : Algebraic Coding Theory 

c 

CBOE GIG GGIOR GO GIGI GIGS GIGI OCG OIOGIGIG IID GIO OIG I I CK 
c 

subroutine dectobin(vectorsize, decinteger, binaryvector) 

c 

integer*4 vectorsize, decinteger, binaryvector(vectorsize) 

integer*4 i, temp1, twostable(6) 

c 

data twostable/1,2,4,8,16,32/ 

c 


c Check to see if the decimal integer is larger than the routine can 
c convert. This can be easily extended by adding to the twostable. 
c 


if (decinteger .gt. 63) 
& print *,’dectobin: decinteger too large’, decinteger 
templ=vectorsize 
do 10 i=1,vectorsize 
if (decinteger.ge.twostable(temp1))then 
binaryvector(temp1)=1 
decinteger=decinteger-twostable(temp1) 
else 
binaryvector(temp1)=0 
end if 
temp1=temp1-1 
10 continue 
return 


c ROUTINE 





annnannanananaanoanaonoeaqanoaenoaaaqana 


anaoaaqanoanaonanaoaqaenanoanqnoo0oaonoendnoaqaenqneaaaqaoeaanga 





matrixgen 
FUNCTION 


This routine is used to generate the H matrix and 

syndrome table necessary for Hamming encode and decode. This 
routine should be called once before calling encodeham and 
decodehan. 


SYNOPSIS 
subroutine matrixgen(codelengthi,codelength2, 
hmatrix,syndrometable) 


formal 


data 1/0 
name type type function 
codelengthi int i number of data bits (63) 
codelength2 int i number of information bits (57) 
hmatrix int o vector to encode an decode by 
syndrometable int o table containing error masks 


DESCRIPTION 


This subroutine is part of a set of subroutines which perform 

a Generalized Hamming Code. As you know, Hamming codes are perfect 
codes and can only detect and correct one error. We added an overall 
parity checkbit, which allows us to detect 2 errors. When 2 errors 
are detected, (in subroutine decodeham.f) no correction attempt is 
made. This would most likely result in more errors. Instead, a flag 
is sent to the calling program notifying it of multiple errors so 

that smoothing may be attempted. The Hamming codes presently supported 
by the routines are (63,57), (31,26), (15,11), and shortened variations 
thereof. It could be made even more general by making minor modifications 
to the dectobin.f subroutine. This routine at present will calculate 

a maximum of 6 bits. 


Hamming routines consist of the following files: 
matrixgen - generates the hmatrix and sydrometable. 


dectobin - does a simple decimal to binary conversion. 
encodeham - generates the codeword and overall paritybit. 
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c decodeham - recovers infobits, checks for errors, corrects 1 

c error, and sends out flag for smoothing. 

c 

c This routine is initializes all of the tables necessary to perform 

c the Hamming code (G Matrix, Syndrome Table) . 
c 

CHS SS SS SS SSS SS SSS SS SSS SS SSS SS SS SSS SSS SSSSSS SSS SSS S SSS SSSSS SSS SSS SSSSSsesz=== 
¢c 

c REFERENCES 

c 

c Lin and Costello : Error Control Coding 

c Berlekamp : Algebraic Coding Theory 

c 


CREO Oa EO oo a naa pana ia a OSS GR aa IG Oa IO ani IK 
c 
subroutine matrixgen(codelength1,codelength2,hmatrix,syndrometable) 
c 
integer*4 codelength1,codelength2 
integer*4 hmatrix(codelength1) ,syndrometable(codelength1i) 
integer*4 itemplate(6) ,ptemplate(57),i,temp1 


This is the data necessary to construct the G Matrix and the Syndrome 
Table. If a larger code is desired, this table can be easily added to. 
All other routines, except the syndrome table construction, 

are general enough to calculate any size Hamming Code. 


aaaqaaaa 


data itemplate/1,2,4,8,16,32/ 

data ptemplate/3,5,6,7,9,10,11,12,13,14,15,17,18,19, 
+ 20,21,22,23,24,25,26,27,28,29,30,31,33,34,35,36,37,38,39,40,41, : 
+ 42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62, 


+ 63/ 
c 
c¢ Construct the parity portion of the hmatrix 
c 


do 30 i=1,codelength2 
hmatrix(i)=ptemplate(i) 

30 continue 

c 

c Construct the identity portion of the hmatrix. 

c 

do 20 i=1, (codelengthi-codelength2) 
hmatrix((codelength2+i))=itemplate(i) 

20 continue : 

c 

c Construct the syndrometable. This routine is rather simple because 
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aqagngnaaaa 








I chose to arrange my G matrix sequentially (Berlrkamp method). 
I placed the parity bits in front in ascending order then added the 


bits left over in ascending order. 


Since our code is linear I can get 


away with this. If a larger Hamming code is needed, then a new 
exception must be generated for each parity bit. 


temp1=1 
do 10 i=1,codelength1 


c 
10 continue 
c 
return 
end 


if(i.eq.1)then 
syndrometable(i)=codelength2+1 
goto 10 

end if 

if(i.eq.2)then 
syndrometable(i)=codelength2+2 
goto 10 

end if 

if (i.eq.4)then 
syndrometable(i)=codelength2+3 
goto 10 

end if 

if (i.eq.8)then 
syndrometable(i)=codelength2+4 
goto 10 

end if 

if (i.eq.16)then 
syndrometable(i)=codelength2+5 
goto 10 

end if 

if (i.eq.32)then 
syndrometable(i)=codelength2+6 
goto 10 

end if 

syndrometable(i)=temp1 

tempi=tempi+t1 





C. DATA INPUT 


Jaaodacooaanoeoee: Gaussian Input eeeedocuoodociiccs 


20. 60. 100. 140. 181. 221. 261. 301. 
341. 381. 421. 462. 502. 542. §82. 622. 
662. 702. 743. 783. 824. 864. 904. 944. 
985. 1025. 1065. 1106. 1146. 1187. 1227. 1268. 

1309. 1350. 1390. 1431. 1472. 1513. 1553. 1594. 
1635. 1676. 1717. 1758. 1798. 1840. 1881. 1922. 
1963. 2005. 2046. 2088. 2129. 2170. 2212. 2254. 
2296. 2337. 2379. 2421. 2463. 2505. 2546. 2589. 
2631. 2673. 2716. 2758. 2800. 2843. 2886. 2928. 
2972. 3014. 3057. 3100. 3144. 3186. 3230. 3273. 
3316. 3360. 3404. 3448. 3492. 3535. 3579. 3624. 
3668. 3712. 3757. 3801. 3846. 3891. 3936. 3981. 
4026. 4072. 4117. 4163. 4209. 4254. 4300. 4346. 
4392. 4439. 4485. 4532. 4579. 4626. 4673. 4720. 
4767. 4815. 4863. 4911. 4959. 5007. 5056. 5104. 
§153. 5201. 5251. 5301. 5350. 5400. 5450. 5500. 
5550. 5600. 5651. §702. 5754. 5805. 5857. 5908. 
5961. 6013. 6066. 6118. 6171. 6225. 6279. 6333. 
6387. 6442. 6496. 6551. 6606. 6662. 6719. 6774. 
6831. 6888. 6945. 7003. 7061. 7119. 7178. 7237. 
7296. 7356. 7417. 7477. 7538. 7599. 7661. 7724. 
7787. 7850. 7914. 7978. 8042. 8108. 8173. 8240. 
8306. 8373. 8441. 8510. 8579. 8649. 8719. 8790. 
8862. 8935. 9007. 9081. 9156. 9231. 9306. 9383. 
9461. 9540. 9619. 9700. 9781. 9863. 9947. 10031. 
10117. 10203. 10290. 10380. 10470. 10561. 10654. 10748. 
10843. 10941. 11039. 11140. 11241. 11345. 11451. 11558. 
11668. 11780. 11893. 12010. 12129. 12250. 12374. 12500. 
12631. 12764. 12901. 13041. 13186. 13334. 13488. 13646. 
13809. 13978. 14153. 14335. 14524. 14721. 14928. 15143. 
15371. 15610. 15864. 16134. 16423. 16734. 17071. 17439. 
17847. 18305. 18827. 19440. 20186. 21150. 22545. 25295. 
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11 
11592 
11559 
20014 
11526 
21965 

756 
16871 
11493 
23980 
23086 

675 
16836 
11588 
19845 
22918 
11460 
21101 
22922 
22920 
19849 
22113 
19847 
23842 
16840 
16842 
11555 
21701 
18756 
23183 
21829 
17762 
11427 
21798 
24100 
21734 
21833 
12609 
21831 
19938 
18760 


362 
11560 
16869 

353 
17932 
24200 
24202 
11525 
16873 
23105 
21763 
12385 
23147 
11556 
24065 
18755 
23714 
18978 

716 
11428 
23087 
24099 
20834 
16839 
23778 
16841 
24012 
11491 

580 
16804 

677 
23085 
20013 
21830 
16867 
12513 

681 
12577 

679 
18759 

584 


Input 


SIO IO OR gk io Hk 


11893 
359 
18563 
11558 
21667 
24201 
11589 
19844 
§27 
11461 
21795 
19874 
21731 
11524 
18691 
21997 
22919 
20015 
18627 
19848 
19850 
16838 
738 
19846 
22917 
355 
18660 
23219 
18724 
19906 
24176 
24195 
21766 
12481 
18596 
21832 
21 
12545 
18695 
23182 
18728 
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360 
11591 
24198 
20033 
11494 
15425 
11557 
21699 
22050 
23682 
21827 
16870 

356 
18659 
14401 
18723 
23201 
22921 
24196 
18893 
23181 
18595 
20961 
21635 

526 
11587 
17889 
20897 
11459 
21797 
18926 
21733 
24272 
22978 
14763 

680 

385 
21799 
16806 
18727 

586 


328 
295 
21889 
262 
10470 
24203 
10533 
229 
23142 
9413 
21859 
23146 
324 
9476 
24173 
196 
15817 
22085 
15815 
8356 
17868 
20802 
20016 
22089 
23140 
291 
21104 
8419 
10435 
23116 
546 
163 
15746 
14760 
14762 
21028 
15780 
14758 
22689 
7299 
587 


18 
$25 
21668 
18758 
21004 
19843 
24106 
23154 
21796 
24102 
21732 
11522 
21100 
16771 
21006 
23120 
481 
21059 
22949 
689 
740 
16737 
19882 
11394 
21964 
449 
20837 
18985 
23781 
18983 
23846 
18924 
20841 
23881 
23848 
16834 
22017 
16775 
17766 
19912 
19914 
17770 
17768 
13441 
718 





18761 
16810 
22146 

§82 
24103 
21700 
18757 
23650 
21828 
21102 

354 
22953 
14465 
22951 
19907 
15617 
24194 
17793 
22056 
22058 
19879 
21602 
22054 
23749 
18982 
14657 
23877 
23813 
14529 
22128 
13409 
23782 
20842 

745 
23879 
23815 
21061 
22947 
16777 
18860 
19913 

555 
21065 
18755 
21063 


18793 
21029 
18781 
609 
705 
19970 
18789 
14722 
21860 
21031 
322 
9474 
17603 
194 
22162 
17731 
17667 
8354 
18849 
20804 
23119 
20974 
23138 
9346 
20709 
13702 
23909 
18987 
17539 
5250 
20778 
20808 
20 
22127 
23911 
20806 
17638 
7266 
22125 
17735 
17671 
24133 
18752 
23153 
17605 


585 
24044 
18694 
19809 
22881 
24105 

581 
16805 

676 
16835 
11586 
22954 
24242 
11458 
18585 
19877 
15649 
22055 
22979 
22057 
23844 
19881 
15521 
14561 
22948 
24275 

741 

§93 
22914 
23878 
23814 
14497 

651 
23880 

743 
23011 
19940 
24045 
18562 
24257 
21666 
19911 
22983 
21064 
22985 


6275 
23042 
19713 
19777 
21030 
18782 
14756 
22721 
15786 
15784 

290 
22955 

8418 
10434 
16545 

547 
162 
15747 
24268 
644 
16673 
22028 
16609 
18956 

7298 
20773 
23187 
15810 

428 
23910 
20807 

648 
23914 
23912 
22082 

646 

8290 
17702 
16779 

551 
17737 
17704 
21067 
18753 
24134 


16809 
18662 
16807 
18726 
13377 
18693 
24104 
21764 
12417 
18629 
11554 

717 
22053 
23716 
18980 
17763 
11426 
20048 
19878 
20836 
23180 
23780 
19880 
19939 
14625 
18530 
18984 
18986 
13707 

742 
19842 
20840 

746 

744 
19876 
20838 
19910 
22981 

§24 
17767 
22148 
15457 

658 
21698 
22051 


22124 
19681 
17538 
§251 
24174 
24107 
20578 
657 
493 
19617 
10530 
226 
24078 
9410 
24172 
13700 
10402 
16577 
20740 
6306 
17869 
19883 
20676 
8322 
23215 
20805 
21921 
13704 
13706 
17901 
647 
20809 
747 
16513 
20742 
98 
15749 
17734 
17670 
19915 
553 
17736 
24077 
19972 
18751 


17857 
16808 
22915 
21998 
18661 
783 
18725 
21035 
15787 
23010 
22952 
11490 
15553 
21996 
22950 
18681 
23748 
24046 
16866 
23876 
23812 
621 
18892 
21007 
17933 
14593 
23845 
22147 
15489 
23216 
20839 
22052 
23 
18981 
23847 
11362 
24209 
16776 
16778 
622 
17769 
16774 
21066 
22986 
22984 


16811 
18811 
19745 
22817 
18959 
22753 
21032 
21034 
22083 
24076 
258 
10466 
20047 
17635 
7362 
18713 
9378 
20708 
22059 
23908 
§13 
17571 
130 
20001 
14689 
645 
13703 
20046 
13705 
20743 
19971 
649 
14723 
20710 
13701 
10338 
24228 
550 
6242 
23044 
554 
552 
15754 
22987 
14754 





16773 
21762 
12353 
21005 
19946 
19944 
18754 
23715 
18979 
17764 
13569 
13633 
24098 
20835 
13505 
23779 

719 
20038 
11489 
15554 
20042 
20040 
15682 
20936 
23111 
16865 
12515 
20934 
12579 
17828 
23206 
22117 
14594 
20036 
24141 
15490 
20966 
24069 
24043 
23617 
20968 
17832 
11361 
23148 
20906 


20611 
24135 
17701 
24137 
19947 
19973 
18786 
14725 
23049 
13699 
22097 
24229 
20739 
14729 
23045 
14727 
20675 
257 
10465 
21898 
20043 
7361 
15714 
9377 
19949 
24036 
22063 
514 
750 
129 
20002 
14690 
21892 
23211 
656 
611 
707 
24040 
24042 
20971 
18825 
24038 
10337 
17899 
20907 


24161 
19909 
22982 
23012 

623 
22150 

578 
16802 
20865 
22152 
23747 
13601 
21932 
23875 
23811 
23016 
23018 
11521 
24273 
23109 

626 
124581 

482 
24237 
20937 
23113 
23053 
16738 
24011 
11393 
24073 

450 
24071 

725 
21003 
19811 
22883 
22122 
22120 

655 
17831 
16833 
22018 

559 

653 


23213 
15750 
17733 
17669 
24233 
23937 
24231 
17636 
22155 
23047 
14726 
20707 
19976 
23907 
14728 
17572 
1 
10497 
9441 
684 
8385 
20996 
19587 
22659 
7329 
17923 
16642 
16706 
22157 
10369 
21069 
18991 
6273 
18821 
19715 
19779 
24039 
22123 
24041 
22723 
23074 
18823 
9313 
18882 
23170 


17935 
21794 
19875 
21730 
13473 
18690 
16772 
19942 

692 
18626 
21060 
17934 
13537 

739 
23015 
18945 

353 
20933 
14466 
21899 
23204 
18618 
24193 
17794 
12483 
20935 
20964 
21601 
12547 
15788 
24074 
14658 
23208 
23210 
14530 
17829 
13410 

691 
23108 
20970 
12419 
17833 
17896 
17898 
24238 


21089 
16481 
549 
4194 
23046 
23117 
20643 
23050 
17700 
19521 
22593 
624 
20771 
460 
20045 
24132 
321 
9473 
21896 
193 
23971 
23021 
21894 
8353 
18850 
23115 
24007 
12643 
23137 
9345 
24075 
24037 
20999 
19683 
21001 
5249 
18822 
16812 
20545 
18826 
20997 
19619 
22692 
7265 
17925 


23683 
21826 
17765 
23084 
18658 
14433 
18722 
19908 
22154 
23013 
19941 
19979 
18594 
14731 
21634 
23017 
11585 
17891 
20899 
11457 
15586 
22064 
15650 
20938 
22116 

722 
20037 

387 
15522 
14562 

752 
21036 
17859 
23209 
22913 
18992 
20932 
14498 
17830 
23205 
22118 
17834 
20903 
12485 
18561 


14724 
21858 
754 
20941 
757 
17604 
18913 
720 
17732 
17668 
13665 
19978 
20803 
14730 
19974 
21953 
289 
21895 
8417 
10433 
16546 
24005 
161 
20939 
21105 
21893 
16674 
24009 
16610 
22691 
7297 
22989 
24108 
15809 
19747 
22819 
22062 
22755 
19917 
588 
22081 
17835 
8289 
19950 
19652 





21062 

674 
24139 
19945 
24235 
19943 
22151 
23051 
22153 
24013 

417 
22980 
23014 
23843 
22149 
20012 
11553 
24068 
19981 
20041 
14764 
20039 
11425 

654 
23112 
23114 
24177 
12611 
23110 
24072 
14626 
18629 
20965 
24070 
23207 
22119 
19841 
22121 
20969 
18827 
20967 

557 
17897 
20905 
17895 


17637 
24136 
24138 
24232 
24234 
17902 
23048 
24230 
548 
19975 
15748 
19977 
643 
21026 
15778 
23019 
10529 
225 
21897 
9409 
23185 
18990 
10401 
16578 
24006 
6305 
18820 
24010 
24008 
8321 
20998 
19651 
21922 
21002 
21000 
22851 
22787 
751 
18824 
16514 
24004 
97 
23974 
24144 
21891 








20904 
21072 
12613 
19812 
22884 

390 
17931 
12583 

673 

392 
17862 
17892 
20900 
17866 
17864 
19816 
19818 

418 
22115 
22160 
22886 

388 
24276 
14564 
22945 
15686 
15623 
16741 
15656 
23181 
20930 
14500 
24163 
16745 
18621 
16743 
14663 

457 
15556 
14566 
20867 
18684 
23745 
13603 
14661 


6241 
749 
19716 
19780 
22130 
17928 
14753 
22724 
395 
17926 
22853 
22789 
19750 
17867 
615 
19784 
22857 
22824 
22855 
22791 
21025 
15777 
613 
24207 
7233 
15718 
16645 
16709 
15722 
18855 
518 
18857 
16680 
16713 
16647 
16711 
20010 
23939 
14696 
685 
23077 
15716 
24178 
20705 
24034 


12517 
23979 
12581 
23054 
13442 
12614 
23107 
21761 
12 
12616 
22092 
24067 
18753 
558 
19815 
24239 
13570 
13634 
24097 
15820 
13506 
12612 
19813 
17797 
14628 
486 
24259 
15655 
430 
17801 
16742 
17799 
22094 
15589 
19873 
15653 
13475 
14665 
456 
21931 
14630 
484 
14565 
22095 
13539 


23976 
19684 
17473 
5217 
17927 
21986 
17929 
23022 
12650 
19620 
709 
22061 
18785 
19719 
19783 
617 
713 
22856 
711 
24240 
18819 
19717 
19781 
8257 
18853 
24109 
21924 
24243 
491 
16646 
16710 
19916 
21091 
16 
24002 
4161 
625 
459 
20641 
24206 
20994 
19554 
22626 
23078 
20769 


17860 
17894 
20902 
12518 
20931 
12582 
24162 
12581 
394 
12584 
22885 
23055 
877 
16801 
20866 
19817 
22889 
13602 
22887 
12516 
19983 
12580 
17827 
24140 
15622 
14596 
20034 
15688 
15 
17802 
17800 
22049 
23681 
523 
23220 
14631 
20011 
14 
14664 
23202 
18620 
14598 
17796 
14629 
23052 


22158 
23977 
19748 
22820 
24204 
22756 
12647 
23973 
12649 
21071 
16844 
23938 

614 
22725 
19752 
19785 

714 

712 
24035 
22822 


24258 
21665 
12549 
24143 
17893 
20901 

391 
21793 

393 
21729 
13474 
18689 
19814 
17865 
12452 
18625 
22890 
22888 
13538 
17861 


19685, 24142 


22758 
65 
20004 
14692 
21890 
16677 
15720 
16613 
17803 
18858 
18856 
520 
522 
16679 
20008 
23020 
14698 
18915 
20006 
14694 
24274 
13667 
18852 
20801 


18946 
11329 
15590 

452 
15654 
15687 

488 
15689 

556 
24208 
15557 
16744 
16746 
15685 
23186 
14599 

458 
19980 
15588 

454 
15652 

419 
22114 
14597 
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16876 
12645 
23975 
22852 
22788 
17930 
21090 
16482 
24003 
4193 
22159 
19718 
19782 
19751 
20995 
19 

22 
21070 
17924 
22854 
22790 
19749 
10305 
16581 
23152 
6209 
15719 
21068 
18721 
16678 
589 
16614 
16712 
23076 
15717 
14695 
20009 
22096 
21929 
16548 
21927 
22156 
20005 
14693 
641 


20035 

389 
15458 
19982 
21697 
12550 
12615 
21825 
12617 
24211 
18657 
14434 
18721 
23203 
17863 

619 

715 
12484 
18593 
20963 
21633 
12548 
22019 
22163 
14660 
21037 

487 
14532 

489 
13443 
18859 
23106 
17798 
12386 

485 

455 
24066 
14632 
23713 
18977 
14662 
13571 
13635 

453 
20833 


23978 
23218 
612 
708 
12646 
560 
18989 
21857 
12648 
22821 
19686 
22757 
18914 
23972 
616 
618 
13666 
22823 
19653 
710 
12644 
21954 
9281 
18883 
23171 
22988 
517 
17505 
5185 
24277 
21987 
20609 
$21 
18854 
519 
21926 
14697 
20007 
21930 
21928 
13697 
16878 
16580 
20737 
23079 


23809 
14533 
17826 
13476 
22023 
16769 
24266 
12450 

483 
21057 
24167 
13540 
21099 
16739 
18948 
13573 
13637 

451 
23217 
13509 
23149 
19810 
22882 
13575 
13639 
22129 
20869 
13542 


516 
19948 
3137 
18889 
21989 
18887 
22060 
17697 
19586 
22658 
21993 
17922 
16643 
16707 
24129 
23943 
24225 
23945 
591 
23041 
19714 
19778 
427 
23941 
18988 
22722 
23075 
686 


16740 
18947 
22026 
22024 
14467 
24165 
19905 
15619 
22093 
17795 
12482 
24169 
24261 
21569 
12546 
22020 
13605 
14659 
18923 
20873 
14531 
20871 

13 
21039 
13607 
18951 
12418 
21963 


16708 
23080 
22027 
18890 
17601 
18916 
23970 
17729 
17665 
13668 
18851 
21992 
21097 
12642 
21956 
21040 
18884 
23172 
18922 
19682 
17537 
13672 
13674 
21988 
20577 
13670 
21960 
19618 


20962 
23150 
21106 

688 
17890 
20898 
24265 
15587 
24263 
15651 

420 
22977 
24170 
15821 

386 
15523 
14563 
20870 
13541 
17858 
20874 
18949 

424 

426 
14499 
24164 

422 
18953 
23009 





21925 

721 
21955 
19884 

590 

724 
21093 
16547 

545 
24110 
15745 

687 
24171 
16675 
21095 
16611 
22690 
23946 
23944 
18921 
20875 
19746 
22818 
13673 
22754 
18917 
21959 
242058 
21961 


23841 
23083 
15524 
24262 
18891 
23179 
22022 
24264 
17761 
13572 
13636 
21995 

755 
13508 
12610 
21038 
19937 
14627 
18497 
24260 

652 
24112 
13640 

425 
18950 
13574 
13638 
18954 
18952 
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23082 
16612 
18888 
23176 
23178 
24267 
18886 
23174 
21094 
16579 
21994 
21098 
21096 
21990 
22956 
13669 

748 
19650 
21923 
18919 
21957 
22850 
22786 
13671 
21092 
16515 
18955 


23777 
24236 
22025 
20049 
15555 

758 
20868 
15683 
22021 
13604 
24168 
12514 
23121 
12578 
24166 

421 
23947 
14595 
20872 
22145 
15491 
13608 
13641 

423 
23649 
13606 

592 





23081 
20673 
23175 
23940 
23177 
17633 
16845 
15715 
21991 
18885 
23173 

659 

§15 
17569 
18918 
20003 
14691 
16877 
18920 
23942 

610 

706 
19969 
21958 
14721 
24111 
21962 
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APPENDIX B: Simulation Data 


INMARSAT Channel 


(AO OOo ORC adagoiiokiioiaiiainiok 


# errors corrected = 2 
there were no errors in QRC cw 


# errors corrected = 2 
success = 763 

you just completed superframe # 999 
# errors corrected = 2 

# errors corrected = 4 
fail = 236 
total errors in sf = 8 

number of 10ms fades in 1000 runs = 1842 
number of 20ms fades in 1000 runs = 254 
number of 40ms fades in 1000 runs = 103 
number of 100ms fades in 1000 runs = 87 
number of 200ms fades in 1000 runs = 24 
you just completed superframe # 1000 
jallerr= 3955 itotbitct= 143856 n= 1000 
the total number of failures = 236 

the total number of successes = 763 

This was a nogolay run 

This was a QRC BE count run 

This run used fading (1=yes,O=no) 1 

s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 99.0000 


ekekkkkekeeeekee Soft Decision Golay **+ ek 


success = 746 
you just completed superframe # 


success = 747 

number of 10ms fades in 1000 
number of 20ms fades in 1000 
number of 40ms fades in 1000 
number of 100ms fades in 1000 
number of 200ms fades in 1000 


you just completed superframe # 
the total number of failures = 


999 
runs = 1909 
runs = 239 
runs = 129 
runs 86 
runs = 32 
1000 
252 


the total number of successes = 747 
This was a soft golay run 

This run used fading (1=yes,0=no) 1 
s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 99.0000 
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B. Gaussian Channel 


so OK NO Code ER KRERKERKKEKK 


errors in 100 bits = 18 

errors in 3000 bits = 337 
errors in 3000 bits = 352 
errors in 3000 bits = 364 
errors in 3000 bits = 367 


you just completed superframe # 200 
the total number of failures = 198 
the total number of successes = 0 
This was a nogolay run 

This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = -24.0000 


s/n during non-fade (dB) = 3.00000 
errors in 100 bits = 11 
errors in 3000 bits = 249 
errors in 3000 bits = 268 
errors in 3000 bits = 265 


errors in 3000 bits = 262 

you just completed superframe # 200 
the total number of failures = 196 
the total number of successes = 1 
This was a nogolay run 

This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = -24.0000 


s/n during non-fade (dB) = 4.00000 
errors in 100 bits = 7 

errors in 3000 bits = 179 

errors in 3000 bits = 182 

errors in 3000 bits = 175 

errors in 3000 bits = 176 

you just completed superframe # 200 
the total number of failures = 197 
the total number of successes = 2 


This was a nogolay run 

This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 5.00000 
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errors in 100 bits = 4 

errors in 3000 bits = 125 

errors in 3000 bits = 128 

errors in 3000 bits = 103 

errors in 3000 bits = 116 

you just completed superframe # 200 
the total number of failures = 192 
the total number of successes = 7 


This was a nogolay run 


This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = 


s/n during non-fade (dB) = 


-24.0000 
6.00000 
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eeeeeeeee Hard Decision Golay ****##k#* 


errors in 100 bits = 14 
errors in 3000 bits = 
errors in 3000 bits = 277 

errors in 3000 bits = 268 

errors in 3000 bits = 278 

you just completed superframe # 200 
the total number of failures = 171 
the total number of successes = 26 
This was a hard golay run 

This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = ~-24.0000 

s/n during non-fade (dB) = 3.00000 


U 
N 
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errors in 100 bits = 1 


errors in 3000 bits = 111 
errors in 3000 bits = 146 
errors in 3000 bits = 145 


errors in 3000 bits = 136 

you just completed superframe # 200 
the total number of failures = 123 
the total number of successes = 72 
This was a hard golay run 

This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = ~-24.0000 

s/n during non-fade (dB) = 4.00000 


errors in 100 bits = 0 

errors in 3000 bits = 41 
errors in 3000 bits = 52 
errors in 3000 bits = 48 


errors in 3000 bits = 36 

you just completed superframe # 200 
the total number of failures = 57 
the total number of successes = 140 
This was a hard golay run 

This run used fading (i=yes,0=no) 0 
s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 5.00000 
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errors in 

errors in 3000 5 

errors in 3000 i 0 

errors in 3000 4 

errors in 3000 i 7 

you just completed superframe # 200 
the total number of failures = 13 


the total number of successes = 186 
This was a hard golay run 
This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 6.00000 





eeekeeeeee Soft Decision Golay ******* 


errors 
errors 
errors 
errors 
errors 


you just completed superframe # 
the total number of failures = 
the total number of successes = 
This was a soft golay run 


in 
in 
in 
in 
in 


100 

3000 
3000 
3000 
3000 


bits 
bits 
bits 
bits 
bits 


10 
146 
192 
172 
182 
200 
135 
62 


This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = 


errors 
errors 
errors 
errors 
errors 


you just completed superframe # 
the total number of failures = 
the total number of successes = 
This was a soft golay run 


in 
in 
in 
in 
in 


100 

3000 
3000 
3000 
3000 


bits 
bits 


bits = 
bits = 


bits 


-24.0000 
s/n during non-fade (dB) = 


3.00000 


68 
200 
77 


This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = 


errors 
errors 
errors 
errors 
errors 


you just completed superframe # 
the total number of failures = 
the total number of successes = 
This was a soft golay run 


in 
in 
in 
in 
in 


100 

3000 
3000 
3000 
3000 


bits 
bits 
bits 
bits 
bits 


-24.0000 
s/n during non-fade (dB) = 


4.00000 


10 
200 
24 


This run used fading (1*yes,0=no) 0 
s/n during fade (dB) = 


-24.0000 
s/n during non-fade (dB) = 


5.00000 
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errors in 100 bits = QO 
errors in 3000 bits = 
etrors in 3000 bits 
errors in 3000 bits = 
errors in 3000 bits = 7 

you just completed superframe # 200 
the total number of failures = 5 
the total number of successes = 194 
This was a soft golay run 

This run used fading (1=yes,0=no) 0 
s/n during fade (dB) = ~-24.0000 

s/n during non-fade (dB) = 6.00000 


tt 
ooo 
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Cc. Constant Burst Width Channel 


KEY: 


fade widths (120, 240, 480, 1200, 2400) correspond to 
burst widths (10, 20, 40, 100, 200) milliseconds. 





cE REE NO Code * kK KK KR RE 
errors in 100 =bits = 4 
errors in 3000 bits = 101 
errors in 3000 bits = 95 
errors in 3000 bits = 86 
errors in 3000 bits = 70 
you just completed superframe # 200 
the total number of failures = 106 
the total number of successes = 93 
This was a nogolay run 
This run used fading (1isyes,0=no) 1 


fading width, bwidth = 120 
s/n during fade (dB) = -24.0000 


s/n during non-fade (dB) = 99.0000 


errors in 100 bits = 7 

errors in 3000 bits = 178 
errors in 3000 bits = 190 
errors in 3000 bits = 188 
errors in 3000 bits = 102 


you just completed superframe # 200 
the total number of failures = 114 
the total number of successes = 85 
This was a nogolay run 

This run used fading (1=yes,0=no) 1 
fading width, bwidth = 240 

s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 0 

errors in 3000 bits = 283 

errors in 3000 bits = 265 

errors in 3000 bits = 324 

errors in 3000 bits = 218 

you just completed superframe # 200 
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the total number of failures = 123 
the total number of successes = 74 
This was a nogolay run 

This run used fading (1=yes,0=no) 1 
fading width, bwidth = 480 

s/n during fade (dB) = -24.0000 


s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 29 
errors in 3000 bits = 472 
errors in 3000 bits = 606 
errors in 3000 =bits = 482 
errors in 3000 bits = 462 
you just completed superframe # 200 
the total number of failures = 135 
the total number of successes = 62 
This was a nogolay run 
This run used fading (1=yes,0=no) 1 
fading width, bwidth = 1200 
s/n during fade (dB) = -24.0000 
s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 26 
errors in 3000 bits = 798 
errors in 3000 bits = 816 
errors in 3000 bits = 852 
errors in 3000 bits = 883 
you just completed superframe # 200 
the total number of failures = 156 
the total number of successes = 39 


This was a nogolay run 
This run used fading (1=yes,0=no) 1 


fading width, bwidth = 2400 
s/n during fade (dB) = -24.0000 
s/n during non-fade (dB) = 99.0000 


eeeaeeeeeeee Hard Decision Golay *****#tee% 


errors in 100 bits = 0 
errors in 3000 bits = 10 
errors in 3000 bits = 13 
errors in 3000 bits = 7 


errors in 3000 bits = 13 

you just completed superframe # 200 
the total number of failures = 11 
the total number of successes = 188 
This was a hard golay run 

This run used fading (i=yes,0=no) 1 


fading width, bwidth = 120 
s/n during fade (dB) = -24.0000 
s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 0 


uN 
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errors in 3000 bits 
errors in 3000 bits = 
errors in 3000 bits 
errors in 3000 bits 18 

you just completed superframe # 200 
the total number of failures = 41 
the total number of successes = 158 
This was a hard golay run 

This run used fading (1=yes,0=no) 1 


“on 
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fading width, bwidth = 240 
s/n during fade (dB) = -24.0000 
s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 0 


errors in 3000 bits = 270 

errors in 3000. ~=s bits = 257 
errors in 3000 bits = 234 

errors in 3000 bits = 305 

you just completed superframe # 200 
the total number of failures = 102 
the total number of successes = 97 
This was a hard golay run 

This run used fading (1*yes,0=no) 1 
fading width, bwidth = 480 

s/n during fade (dB) = -24.0000 
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s/n during non-fade (dB) = 99.0000 


errors in 100 bits = 0 

errors in 3000 bits = 698 

errors in 3000 bits = 543 

errors in 3000 bits = 733 

errors in 3000 bits = 616 

you just completed superframe # 200 
the total number of failures = 151 
the total number of successes = 48 


This was a hard golay run 
This run used fading (1=yes,0=no) 1 


fading width, bwidth = 1200 
s/n during fade (dB) = -24.0000 
s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 14 
errors in 3000 ~=bits = 739 
errors in 3000 bits = 972 


errors in 3000 bits = 986 

you just completed superframe # 200 
the total number of failures = 133 
the total number of successes = 19 
This was a hard golay run 

This run used fading (1=yes,0=no) 1 
fading width, bwidth = 2400 

s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 99.0000 
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eeeeeeeeEEE Soft Decision Golay ******ke* 


errors in 100 bits = 0 
errors in 3000 bits = 0 
errors in 3000 bits = 0 
errors in 3000 bits = 0 
errors in 3000 bits = (¢) 


you just completed superframe # 200 
the total number of failures = 2 
the total number of successes = 197 
This was a soft golay run 

This run used fading (1i=yes,0=no) 1 


fading width, bwidth = 120 

s/n during fade (dB) = -24.0000 
s/n during non-fade (dB) = 99.0000 
errors in 100 bits= 0 

errors in 3000 bits = (0) 

errors in 3000 bits = 76 

errors in 3000 ~=bits = 10 

errors in 3000 bits = 26 


you just completed superframe # 200 
the total number of failures = 21 
the total number of successes = 178 
This was a soft golay run 

This run used fading (i=yes,0=no) 1 
fading width, bwidth = 240 

s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 99.0000 


errors in 100 bits = 0 


errors in 3000 bits = 111 
errors in 3000 bits = 137 
errors in 3000 bits = 133 
errors in 3000 bits = 153 


you just completed superframe # 200 
the total number of failures = 60 
the total number of successes = 139 
This was a soft golay run 

This run used fading (1=yes,0=no) 1 
fading width, bwidth = 480 

s/n during fade (dB) = -24.0000 
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s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 

errors in 3000 bits = 412 

errors in 3000 bits = §21 
errors in 3000 bits = 647 

errors in 3000 bits = 509 

you just completed superframe # 200 
the total number of failures = 133 
the total number of successes = 66 


This was a soft golay run 
This run used fading (1=yes,0=no) 1 


fading width, bwidth = 1200 
s/n during fade (dB) = -24.0000 
s/n during non-fade (dB) = 99.0000 
errors in 100 bits = 
errors in 3000 bits = 769 
errors in 3000 bits = 711 
errors in 3000 bits = 976 
errors in 3000 bits = 757 


you just completed superframe # 200 
the total number of failures = 160 

the total number of successes = 39 

This was a soft golay run 

This run used fading (i=yes,0=no) 1 
fading width, bwidth = 2400 

s/n during fade (dB) = -24.0000 

s/n during non-fade (dB) = 99.0000 








Hamming Codeword Interleaving Table : 


CW4 


731 751 

1001 1021 
1271 1291 
1541 1561 
1811 1831 
2081 2101 
2351 2371 


2621 2641 


CW10 CW1l CW12 CW13 CWl4 CW15 CW16 CW17 CW18 


1151 


1411 1421 


1681 1691 


1951 196} 


2221 2231 


2491 2501 


2761 2771 ‘ 
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Re ye ee ee! 





Golay Codeword Interleaving Table 











CW CW2 CW3 CcW4 CWS 


QRC Codeword Interleaving Table 
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